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Preface 



These notes were used in a short graduate course on branching processes the author gave 
in Beijing Normal University. The following main topics are covered: scaling limits 
of Galton-Watson processes, continuous-state branching processes, extinction probabil- 
ities, conditional limit theorems, decompositions of sample paths, martingale problems, 
stochastic equations, Lamperti's transformations, independent and dependent immigra- 
tion processes. Some of the results are simplified versions of those in the author's book 
"Measure-valued branching Markov processes" (Springer, 2011). We hope these simpli- 
fied results will set out the main ideas in an easy way and lead the reader to a quick access 
of the subject. 



Zenghu Li 
Beijing, China 
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Chapter 1 
Preliminaries 



In this chapter, we discuss the basic properties of Laplace transforms of finite measures 
on the positive half line. In particular, we give some characterizations of the weak con- 
vergence of those measures in terms of their Laplace transforms. Based on these results, 
a general representation for infinitely divisible distributions on the positive half line is 
established. We also give some characterizations of continuous functions on the positive 
half line with Levy-Khintchine type representations. 



1.1 Laplace transforms of measures 

In this section, we discuss the basic properties of Laplace transforms of finite measures 
on the positive half line R+ := [0, oo). Let B(R+) = b=^(M+) be the set of bounded 
Borel functions on M_|_. Given a finite measure (x on M+, we define the Laplace transform 
of /i by 

LA>')= / e-^Xdx), A>0. (1.1.1) 
Jo 

Theorem 1.1.1 A finite measure on IR+ is uniquely determined by its Laplace transform. 

Proof. Suppose that /ii and /i2 are finite measures on M_|_ and /^^^(A) = L^^{X) for all 
A > 0. Let ^ = {x H-^ e"^^ : A > 0} and let ^ be the class of functions F E B{R+) so 
that 

POO POO 

/ F(x)/ii(dx) = / F{x)fi2{dx). 
Jo Jo 

Then is closed under multiplication and ^ is a monotone vector space containing 

J^. It is easy to see cr(J^) = ^(M+). Then the monotone class theorem implies =Sf D 

ha{J(f) = B{R+). That proves the desired result. □ 
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Theorem 1.1.2 Let be finite measures on and let A i— ?■ L{\) be a continuous 
function on [0, oo). If there is a dense subset D o/ (0, oo) to that lim„_^oo -^/i„(A) = 
L(A) for every \ E D, then there is a finite measure fj, on R+ such that L^ = L and 
lim„_^oo fJ'n = fJ' by weak convergence. 

Proof. We can regard each i^in as a finite measure on R+ := [0, oo], tlie one-point com- 
pactification of [0, oo). Let F„ denote the distribution function of By applying Helly's 
theorem one can see that any subsequence of {Fn} contains a weakly convergent subse- 
quence {Fn^^}. Then the corresponding subsequence {/Unfe} converges weakly on M+ to a 
finite measure /i. It follows that 

/i(t+) = lim /i„,(M+) = lim L^„^(0) = L(0). 
Moreover, for A G -D we have 

/•oo 

e-^^/i(dx) = lim / e-^^/i„,(da;) = L(A), (1.1.2) 
Jo 

where e~'^ °° = by convention. By letting A 0+ along D in (11.1.21) and using 
the continuity of L at A = we find /i(M+) = -^^(0), so /i is supported by M+. Then 
lim„_j.oo /infe = weakly on IR+. It is easy to see that (11.1.21) in fact holds for all A > 0, 
so we have L^j = L. By a standard argument one sees lim„_j.oo l^n = 1^ weakly on IR+. □ 



Theorem 1.1.3 Let /ii, yU2, • • • and fi be finite measures on M+. Then /in — > /i weakly if 
and only if L^^^{\) — )■ L^{\) for every A > 0. 

Proof, li ^ IJ^ weakly, we have lim„^oo -^At„(A) = -^^^t(A) for every A > by 
dominated convergence. The converse assertion is a consequence of Theorem 1 1.1.21 □ 

We next give a necessary and sufficient condition for a continuous real function to be 
the Laplace transform of a finite measure on M+. For a constant c > and a function / 
on an interval T C M we write 

AJ(A) = /(A + c)-/(A), A,A + cGT. 

Let A° be the identity and define A" = A"^^ Ac for n>l inductively. Then we have 

i=0 ^ ^ 

The Bernstein polynomials of a function / on [0, 1] are given by 

BfAs) = 5Z P Al/^/(0)s', < s < 1, m = 1, 2, . . . . (1.1.3) 
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It is well-known that Bf^m{s) f{s) uniformly as m — j- oo; see, e.g., Feller (1971, 
p. 222). A real function 6 on [0, oo) is said to be completely monotone if it satisfies 

(-l)*A',e(A) > 0, A>0,c>0,z = 0,l,2,.... (1.1.4) 

Theorem 1.1.4 A continuous real function 9 on [0, oo) is the Laplace transform of a finite 
measure fi on M+ if and only if it is completely monotone. 

Proof. If 9 is the Laplace transform of a finite measure on it is clearly a completely 
monotone function. Conversely, suppose that (11.1.41) holds. For fixed a > 0, we let 
la{s) = 9{a — as) for < s < 1. The complete monotonicity of 9 implies 

Al/^7a(0) >0, z = 0,l,...,m. 

Then the Bernstein polynomial By^ rn{s) has positive coefficients, so B^^^m{e^^^^) is the 
Laplace transform of a finite measure on M_|_. By Theorem ll.l.2[ 



^(A) = lim lim B^^^mie-^^"), A > 0, 

a— )-oo m— )-oo 



is the Laplace transform of a finite measure on M+. □ 

We often use a variation of the Laplace transform in dealing with cr-finite measures on 
(0, oo). A typical case is considered in the following: 

Theorem 1.1.5 Let /ii and /i2 be two a-finite measures on (0, oo). If for every A > 0, 

/•oo /"OO 

/ (l-e-">i(da;)= / (1 - e-"^)/.2(da;) (1.1.5) 
Jo Jo 

and the value is finite, then we have /ii = /i2- 

Proof. By setting Hi{{0}) = /i2({0}) = we extend /ii and /i2 to cr-finite measures on 
[0, oo). Taking the difference of (11.1.51) for A and A + 1 we obtain 

/ e-^"(l-e-")/ii(dx) = / e-^"(l-e-")/i2(dx). 
Jo Jo 

Then the result of Theorem 1 1 . 1 . 1 1 implies that 

(1 - e"^)/ii(dx) = (1 - e"^)/i2(dx) 

as finite measures on [0, oo). Since 1 — e""^ is strictly positive on (0, oo), it follows that 
111 = jj2 as cr-finite measures on (0, oo). □ 
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Now let us consider a complete separable metric space E with the Borel a-algebra 
denoted by 3S{E). Suppose that /i is a strictly positive bounded Borel function on E. Let 
Bh{E) be the set of Borel functions f on E such that |/| < const ■ h. Let be the 
set of Borel measures ji on E such that hdji < oo. Let be the ci-algebra on Mh 
generated by the mappings 

Je 

Given a finite measure Q on {Mh, ^h), we define the Laplace functional Lq of Q by 

LQ{f)= [ e-^(^)g(dz/), feBhiE^. (1.1.6) 

A random element X taking values on {Mh, ^h) is called a random measure on E. The 
Laplace functional of a random measure means the Laplace functional of its distribution 
on {Mh, ^h) ■ The reader may refer to Kallenberg (1975) or Li (201 1) for the basic theory 
of random measure. In particular, the proofs of the following results can be found in the 
two references: 

Theorem 1.1.6 A finite measure on {M^, -Mh) is uniquely determined by its Laplace func- 
tional. 

Suppose that A is a a-finite measure on {E, .!^{E)). A random measure X on is 
called a Poisson random measure with intensity A provided: 

(1) for each B E 3S{E) with A(-B) < oo, the random variable X{B) has the Poisson 
distribution with parameter \{B), that is, 

nX{B)=n} = ^e-'^^\ n = 0,1,2,...; 

(2) if Bi, . . . , Bn E .^{E) are disjoint and \{Bi) < oo for each i = 1, . . . ,n, then 
X{Bi), . . . , X{Bn) are mutually independent random variables. 

Theorem 1.1.7 A random measure X on E is Poissonian with intensity A E Mh{E) if 
and only if its Laplace functional is given by 

Eexp{-X(/)} = exp { - ^(1 - e-^(^))A(dx)}, / E Bh{E)+. (1.1.7) 
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Proof. Suppose that X is a Poisson random measure on E with intensity A. Let Bi, . . . , Bn G 
^(E) be disjoint sets satisfying \{Bi) < oo for each i = 1, . . . ,n. For any constants 
"i, • • • , ttn > we can use the above two properties to see 

n n 

Eexp I - = exp | - ^(1 - e-"0A(5,)}. (1.1.8) 

i=l 1=1 

Then we get (11.1.71) by approximating / G B h{E)^ hy simple functions and using dom- 
inated convergence. Conversely, if the Laplace functional of X is given by (11.1.71) . we 
may apply the equality to the simple function / = XlILi '^i^B^ to get (11.1.81) . Then X 
satisfies the above two properties in the definition of a Poisson random measure on E 
with intensity A. 



1.2 Infinitely divisible distributions 

For probability measures fii and fi2 on M+, the product fii x is a probability measure on 
. The image of fii x fi2 under the mapping (xi ,X2) h-> xi + X2 is called the convolution 
of /ii and /i2 and is denoted by /ii * ^2, which is a probability measure on M+. According 
to the definition, for any F E B(R+) we have 

/•oo /"OO /*oo 

/ F(x)(/ii * /i2)(dx) = / /ii(dxi) / + X2)/i2(dx2). (1-2.1) 

Jo Jo Jo 

Clearly, if and ^2 are independent random variables with distributions /ii and fX2 on M+, 
respectively, then the random variable ^1 + ^2 has distribution /xi * /i2. It is easy to show 
that 

L^,.^,(A) = L^,(A)L^,(A), A > 0. (1.2.2) 

Let = 5o and define /i*" = * inductively for integers n > 1. We say a 

probability distribution /i on IR+ is infinitely divisible if for each integer n > 1, there is a 
probability /i„ such that /i = /i*". In this case, we call //„ the n-th root of /i. A positive 
random variable ^ is said to be infinitely divisible if it has infinitely divisible distribution 
on R+. 

We next give a characterization for the class of infinitely divisible probability mea- 
sures on R_|_. Write ^' E .y if X ^ V'(A) is a positive function on [0, 00) with the 
representation 

POO 

^(A) = hX+ (1 - e-^")/(dM), (1.2.3) 
Jo 

where h > and (1 A u)l{du) is a finite measure on (0, 00). 
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Proposition 1.2.1 The pair {h, I) in (11.2.31) is uniquely determined by the function ip G 
J. 

Proof. Suppose that ^ can also be represented by (11.2.31) with (/i, /) replaced by {h! , /'). 
For A > and 6' > 0, we can evaluate + ^) — '\\){&) with the two representations and 
get 

/•OO /'OO 

/lA + / (1 - e"^")e-^"/(dM) = h'\+ (l - e-^")e~^"/'(dM). 
Jo Jo 

By letting 6* — oo we get h = h', and so l{du) = l'{du) by Theorem 1 1.1. 51 □ 

Theorem 1.2.2 Suppose that ip is a continuous function on [0, oo). If there is a sequence 
{ipn} C such that ipi^X) = lim„^oo "ipnWfor all A > 0, then il) G J^. 

Proof. Suppose that ipn G is given by (11.2.31) with {h,l) replaced by (hnjn)- We 
can define a finite measure F„ on M+ by setting F„({0}) = ^^^({oo}) = and 
Fr,{du) = (1 - e-")/„(dM) for < M < OO. For A > let 

r (1 - e-")-i(l - e-"^) if < M < OO, 

^(n,A) = < A ifw = 0, (1.2.4) 

[^1 if ti = OO. 

Then we have 

V'n(A) = / eKA)F„(dn), A>0. 



It is evident that {F„(]R+)} is abounded sequence. Take any subsequence {Fn^} C {Fn} 
such that limfc^oo Fn^ = F weakly for a finite measure F on M+. Since m h-> ^(m, A) is 
continuous on M+, we have 

^(A) = / ^(n,A)F(du), A > 0. 
Jr+ 

Observe also that Mm.n^ooi'i^/'n) = ?/'(0) = implies F({oo}) = 0. Then the desired 
conclusion follows by a change of the integration variable. □ 



Theorem 1.2.3 The relation ip = — log establishes a one-to-one correspondence be- 
tween the functions ip E ^ and infinitely divisible probability measures fi on M_|_. 
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Proof. Suppose that tp E is, given by (I1.2.3I) . Let be a Poisson random measure on 
(0, oo) with intensity l{du) and let 



^ = h+ xN{dx). 
Jo 

By Theorem ll. 1.71 for any A > we have 

E e-^« = exp I - /lA - ^ (l - e"^") /(du) } . 

Then 'ip = — log L^j for a probability measure yU on M+. Similarly, for each integer n > 1 
there is a probability measure /i„ on M+ so that ip/n = — \ogL^^^. It is easy to see that 
/i™ = ji. That gives the infinite divisibility of ji. Conversely, suppose that ^/^ = — log 
for an infinitely divisible probability measure ji on IR+. For n > 1 let be the ra-th root 
of /i. Then 

poo 

i}{X) = lim n[l - e-""''^^^)] = lim / (l - e-^^)n/i„(da;). 



By Theorem ll ■2.2l we have ip E J^. □ 

The above theorem gives a complete characterization of infinitely divisible probability 
measures on M+. We write fi = I{h, I) if /i is an infinitely divisible probability measure 
on R+ with tp := — log given by (I1.2.3I) . 

Theorem 1.2.4 If^jji, 1^2 E J, then '^\0'^2^ ^■ 

Proof. For every x > we clearly have xi\)2 ^ ^ so there is an infinitely divisible prob- 
ability measure on M+ satisfying — log Ly^ = xil)2. By a monotone class argument one 
can see Uxidy) is a probability kernel on M+. Let ji be the infinitely divisible probability 
measure on M+ with — log = ipi and define 

poo 

V{(^y) = / fJ'{dx)ux(dy), y>0. 
Jo 

It is not hard to show that — logL^ = tpi o tp2. By the same reasoning, for each integer 
n > 1 there is a probability measure ?7„ such that — log = n^^ipi o ?/)2- Then 77 = 77*" 
and hence 77 is infinitely divisible. By Theorem ll .2.31 we conclude that o ^2 ^ D 

Example 1.2.1 Let b > and a > 0. The Gamma distribution 7 on IR+ with parameters 
(6, a) is defined by 
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which reduces to the exponential distribution when 6 = 1. The Laplace transform of 7 is 

It is easily seen that 7 is infinitely divisible and its n-th root is the Gamma distribution 
with parameters {b/n,a). 

Example 1.2.2 For c > and < a < 1 the function A ^ cA" admits the representation 
(|1.2.3I) . Indeed, it is simple to show 

A" = ..77^ - e-^")-^, A > 0. (1.2.5) 



The infinitely divisible probability measure u on M+ satisfying — logLy(A) = cA" is 
known as the one-sided stable distribution with index < a < 1. This distribution does 
not charge zero and is absolutely continuous with respect to the Lebesgue measure on 
(0, 00) with continuous density. For « = 1/2 it has density 

q{x) := ^a;"=^/2e-^'/^", x > 0. 

For a general index the density can be given using an infinite series; see, e.g., Sato (1999, 
p.88). 



1.3 Levy-Khintchine type representations 

In this section, we give some criteria for continuous functions on [0, 00) to have Levy- 
Khintchine type representations. The results are useful in the study of high-density limits 
of discrete branching processes. For u > and A > let 

Uu, A) = e-^" - 1 - (1 + u-)-' J2 n = 1, 2, . . . . 

1=1 

We are interested in functions on [0, 00) with the representation 

n—l 

0(A) = ^a,A^+/ enKA)(l-e-")-"G(dM), A > 0, (1.3.1) 

where n > 1 is an integer, {oq, . . . , a„_i} is a set of constants and G{du) is a finite 
measure on IR+. The value at m = of the integrand in (11.3.11) is defined by continuity as 
(— A)"/n!. The following theorem was proved in Li (1991, 2011): 
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Theorem 1.3.1 A continuous real function on [0, oo) has the representation (11.3.11) if 
and only if for every c > the function 

^,(A) := (-1)"AXA), A>0 (1.3.2) 

is the Laplace transform of a finite measure on M+. 

Based on the above theorem we can give canonical representations for the limit func- 
tions of some sequences involving probability generating functions. Let {a^} be a se- 
quence of positive numbers and let {gk\ be a sequence of probability generating functions, 
that is, 

oo 

9k{z) = ^Vkiz\ \z\ < 1, 

i=0 

where pki > and Xli^o^'^* ^ ^- ^^^^ consider the sequence of functions {ipk} 
defined by 

V'fc(A) = afc[l-(^,(l-A/A;)], < A < A;. (1.3.3) 

Theorem 1.3.2 If the sequence {'ipk\ defined by (|1.3.3I) converges to a continuous real 
function ip on [0, oo), then the limit function belongs to the class defined by (|1.2.3I) . 

Proof For any c, A > and sufficiently large > 1 we have 

A,V^fc(A) = -«fcA,(7fc(l - ■/A;)(A). 

Since for each integer i>\ the z-th derivative (7^ is a power series with positive coeffi- 
cients, we have 

(-l)^^A,V^fc(A) = -k-^a^^J^(\ - ■/k){\) > 0. 

By the mean-value theorem, one sees inductively (— l)*A)jAc^/'fc(A) > 0. Letting — )■ 00 
we obtain (— l)*A^Ac^/'(A) > 0. Then Ac^/'(A) is a completely monotone function of 
A > 0, so by Theorem 1 1.1. 41 it is the Laplace transform of a finite measure on M+. Since 
'ijj{0) = limfc_j.oo ^fc(O) = 0, by Theorem 11.3.11 there is a finite measure F on M_|_ so that 

POO 

^|J{X)= / (l-e"^")(l-e-^)"^F(dn), 
^0 



where the value of the integrand at u = is defined as A by continuity. Then (11.2.31) 
follows with 13 = F{{0}) and n{du) = (1 - e~")-iF(dM) for u>0. □ 
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Example 1.3.1 Suppose that (7 is a probability generating function so that (3 := g'{l~) < 
00. Let ak = k and gkiz) = g{z). Then the sequence ^a;(A) defined by (11.3.31) converges 
to /3A as A; — 00. 

Example 1.3.2 For any < a < 1 the function ^{X) = A" has the representation (11.2.31) . 
For a = 1 that is trivial, and for < a < 1 that follows from (11.2.51) . Let i^ki^) be defined 
by (fT33l) with = A;" and ^fc(^) = 1 - (1 - z)". Then ^/'^(A) = A° for < A < A;. 

In the study of limit theorems of branching models, we shall also need to consider 
the limit of another function sequence defined as follows. Let {ak} and {gk}he given as 
above and let 

0fc(A) = ak[gk{l - \/k) - (1 - A/A;)], < A < A;. (1.3.4) 

Theorem 1.3.3 If the sequence {(j)k} defined by (11.3.41) converges to a continuous real 
function on [0, 00), then the limit function has the representation 

0(A) = aA + cA^ + / (e-^"-l + ^)m(du), (1.3.5) 

where c > and a are constants, and midu) is a a-finite measure on (0, 00) satisfying 

(1 A n^)m(dn) < 00. (1.3.6) 

Proof. Since 0(0) = limfc_^oo 0fc(O) = 0, arguing as in the proof of Theorem 1 1 . 3 . 21 we see 
that has the representation (11.3.11) with n = 2 and = 0, which can be rewritten into 
the equivalent form (11.3.51) . □ 

Now let us consider a special case of the function given by (11.3.51) . Observe that if 
the measure m{du) satisfies the integrability condition 

POD 

/ (n A n^) m(d?i) < 00, (1.3.7) 
Jo 

we have 

POO 

0(A) = 6A + cA^ + / (e-^"-l + AM)m(dM), (1.3.8) 
Jo 

where 



1.3. LEVY-KHINTCHINE TYPE REPRESENTATIONS 



11 



Proposition 1.3.4 A function with the representation (11.3.51) is locally Lipschitz if and 
only ?T (II. 3.71) holds. 

Proof. For computational convenience we first rewrite (|1.3.5I) as 

noo 

</,(A) = 6iA + cA^ + / (e~^"-l + AMl{„<i})m(dM), (1.3.9) 
Jo 

where 




By applying dominated convergence to (11.3.91) . for each A > we have 

0'(A) = hi + 2c\+ I u{l - e~^")m(dM) - [ ue~^''m{du). 
Then we use monotone convergence to the two integrals to get 

0'(O+) = bi- um{du). 

i(l,oo) 

If is locally Lipschitz, we have 0'(O+) > — oo and the integral on the right-hand side 
is finite. This together with (|1.3.6I) implies (11.3.71) . Conversely, if (|1.3.7I) holds, then 0' is 
bounded on each bounded interval and so is locally Lipschitz. □ 



Corollary 1.3.5 If the sequence {0^} defined by (11.3.41) is uniformly Lipschitz on each 
bounded interval and (pki^) (PWfar all X > as k ^ oo, then the limit function has 
the representation (11.3.81) . 



Example 1.3.3 Suppose that g is a probability generating function so that g'{l—) = 1 and 
c := g"{l—)/2 < oo. Let ak = k"^ and gk{z) = g{z). By Taylor's expansion it is easy to 
show that the sequence 0fe(A) defined by (11.3.41) converges to cA^ as A; — t- oo. 



Example 1.3.4 For < « < 1 the function 0(A) = — A° has the representation (11.3.51) . 
That follows from (11.2.51) as we notice 

/•°° / u \ du _ / 1 \du 
io \l + ^)u^ ~ Jo \T+^)^ 

The function is the limit of the sequence 0fc(A) defined by (11.3.41) with ak = k" and 

gkiz) = 1 - (1 - zr. 
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Example 1.3.5 For any 1 < a < 2 the function 0(A) = A" can be represented in the form 
of (I1.3.8I) . In particular, for 1 < a < 2 we have 

A" = ^^^ r(e---l + A.)4^, A>0. 

Let 0fc(A) be defined by (ll.3.4|) with ak = ak"" and gkiz) = z + — 2;)". Then 

0;,,(A) = A" forO < A < A;. 



Chapter 2 

Continuous-state branching processes 



In this chapter, we first give a construction of CB-processes as the scaling limits of discrete 
Galton-Watson branching processes. This approach also gives the interpretations of the 
CB-processes. We shall study some basic properties of the CB-processes. In particular, 
some conditional limit theorems will be given. We also give a reconstruction of the sample 
paths of the CB-processes in terms of excursions. 



2.1 Construction by scaling limits 

Suppose that j : n,i = 1, 2, . . .} is a family of positive integer-valued i.i.d. random 
variables with distribution given by the probability generating function g. Given the pos- 
itive integer a; (0) = m, we define inductively 

x{n—l) 

^in) = J2 n = 1,2,.... (2.1.1) 

i=l 

It is easy to show that {x{n) : n = 0, 1, 2, . . .} is a discrete-time positive integer-valued 
Markov chain with transition matrix P{i,j) defined by 

oo 

Y,P{^,j)z' =9izy, 1 = 0,1,2,..., \z\<l. (2.1.2) 

j=0 

The random variable x{n) can be thought of as the number of individuals in generation 
n > of an evolving particle system. After one unit time, each of the x{n) particles splits 
independently of others into a random number of offspring according to the distribution 
given by g; see, e.g., Athreya and Ney (1972). For n > the n-step transition matrix 
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P'^ii.i) is determined by 

OC' 

5^P"(^,j>^ = (7"(^)\ 2 = 0,1,2,..., 1^1 <1, (2.1.3) 

i=o 

where g'^{z) is defined by g'^{z) = g{g"'~^{z)) successively with g^{z) = z. We call any 
positive integer-valued Markov chain with transition matrix given by (12.1.21) or (12.1.31) a 
Galton-Watson branching process (GW-process). If g'{l—) < oo, the first moment of the 
discrete distribution j); j = 0, 1, 2, . . .} is given by 

oo 

5^JP"(^,J)=^^?'(l-)^ (2.1.4) 

i=i 

which can be obtained by differentiating both sides of (12.1.31) . 

Now suppose we have a sequence of GW-processes {xk{n) : n > 0} with offspring 
distribution given by the sequence of probability generating functions {gk}- Let Zk{;n) = 
Xk{n)/k for n > 0. Then {zk{n) : > 0} is a Markov chain with state space : = 
{0,l/k,2/k, . . .} and n-step transition probability P^ix, dy) determined by 

/ e-^^P,"(x,dy)=^7,"(e-^/Y^ A > 0. (2.1.5) 

Suppose that {7^} is a positive sequence so that 7^. — )■ 00 as A; — 00. Let [7^^] denote the 
integer part of 7fct > 0. We are interested in the asymptotic behavior of the sequence of 
continuous time processes {zk{Ylkt\) : t > 0}. By (12.1.51) we have 

I e^^^Pi^'=*](x,di/) = exp{-xt;,(t,A)}, (2.1.6) 

JEk 

where 

VkitA) = ~k\ogg^^''\e-'/'), A > 0. (2.1.7) 

Clearly, if Zk{0) = x E Ek, then the probability p]^'°*^(x, ■) gives the distribution of 
Zkillkt]) on M^. Let us consider the function sequences 

Gk{z) = k^kigkie-'^") - z > 0, (2.1.8) 

and 

Mz) = k^kMl - z/k) - (1 - z/k)l 0<z<k. (2.1.9) 

Proposition 2.1.1 The sequence {Gk} is uniformly Lipschitz on each bounded interval if 
and only if so is {0fc}. In this case, we have limfc_j.oo — Gk{z)\ = uniformly on 

each bounded interval. 
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Proof. From (12.1.81) and (12.1.91) it is simple to see that 

G',{^) = 7/ce-^/'[l - g',{e~^/% z > 0, (2.1.10) 

and 

(j)i{z) = 7fc[l - g'kil - z/k)], 0<z<k. (2.1.11) 

Clearly, the sequence {G^} is uniformly bounded on each bounded interval if and only 
if so is {(/)'}.}. Then the first assertion is immediate. We next assume {Gk} is uniformly 
Lipschitz on each bounded interval. Let a > 0. By the mean-value theorem, for A; > a 
and < z < a we have 

Gkiz) - Mz) = k7k [gkie-'^") - gkil - z/k) - e-^/^ + (1 - z/k)] 

= kjkWkiVk) - l](e-^/' - 1 + z/k), (2.1.12) 

where 

1 - a/k <1 - z/k <r]k < e""/^' < 1. 
Choose kg > a so that e"^"/'^^ < ^ — a/ko. Then e~^"/^ < 1 — a//;; for /c > /cq and hence 

1k\g'kiVk) - l\ < sup 7fc|^^(e-"/^) - 1|, A;>A;o- 

0<z<2a 

Since {Gk} is uniformly Lipschitz on each bounded interval, the sequence (12.1.101) is 
uniformly bounded on [0, 2a]. Then {'jklgkiVk) — M '■ ^ > kg} is a bounded sequence. 
Now the desired result follows from (12.1.121) . □ 

By the above proposition, if either {Gk} or {(pk} is uniformly Lipschitz on each 
bounded interval, then they converge or diverge simultaneously and in the convergent 
case they have the same limit. For the convenience of statement of the results, we formu- 
late the following conditions: 

Condition 2.1.2 The sequence {Gk} is uniformly Lipschitz on [0, a] for every a > and 
there is a function on [0, oo) so that Gk{z) — (p{z) uniformly on [0, a] for every a > 
as k ^ CO. 

Proposition 2.1.3 Suppose that Condition \2.L2\ is satisfied. Then the function has 
representation 

/•oo 

(f){z) = bz + cz"^ + (e"^" - l + ZM)m(dM), z>0, (2.1.13) 
Jo 

where c > and b are constants and {u A v?)m{(\u) is a finite measure on (0, oo). 
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Proof. By Proposition 12.1 .11 the sequence {(pk} is uniformly Lipschitz on [0,a] and 
<Pkiz) — )■ (f){z) uniformly on [0, a] for every a > 0. Then the result follows by Corol- 
larv [T331 □ 



Proposition 2.1.4 For any function (j) with representation (I2.1.13|) there is a sequence 
{Gk} in the form of (12.1.81) satisfying Condition \2.1.2\ 

Proof. By Proposition l2.1.1l it suffices to construct a sequence {0^} in the form of (12.1.91) 
uniformly Lipschitz on [0, a] and (pkiz) — )■ 0(2;) uniformly on [0, a] for every a > 0. 
To simplify the formulations we decompose the function into two parts. Let (t)o{z) = 
(f){z) — hz. We first define 

/■oo 

7o,fc = (1 + 2c)k + / u{l- e-'=")m(dM) 

and 

QqA'^) = ^ + ^~So^fc0o(/i;(l - 2)), \z\ < 1. 

It is easy to see that z (-)■ go,k{z) is an analytic function in ( — 1, 1) satisfying (^o,fc(l) = 1 
and 

d" 

^5o,fc(0) > 0, n > 0. 

Therefore go,k{ ) is a probability generating function. Let 0o,fe be defined by (12.1.91) with 
ilk, 9k) replaced by (7o,fc, fi'o.fc)- Then 0o,fc(-2) = 0o(-2) for < 2; < A;. That completes the 
proof if 6 = 0. In the case 6 7^ 0, we set 




Let 7i,fe = \b\ and let 01,^(2;) be defined by (12.1.91) with {'jk, Qk) replaced by {-yi^k, 9i,k)- 
Thus we have 

^,^,{z) = bz + ^{\b\-b)z\ 

Finally, let 7fc = 7o,fc + 7i,fc and gk = lk^{lo,kgo,k + 7i,fc5'i,fe)- Then the sequence 0fc(z) 
defined by (12.1.91 ) is equal to (po,k{z) + (pi,k{z) which satisfies the required condition. □ 

Lemma 2.1.5 Suppose that the sequence {Gk} defined by (12.1.81) is uniformly Lipschitz 
on [0, 1]. Then there are constants B,N > such that Vk{t, A) < Xe^^ for every t, \ > 
and k > N. 
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Proof. Let bk := G'f.{0+) for A; > 1. Since {Gk} is uniformly Lipschitz on [0, 1], the 
sequence {bk} is bounded. From (12.1.81) we have b^ = 7fc[l — g[{l—)]. By (12.1.41) it is not 
hard to obtain 

Let i? > be a constant such that 2\bk\ < B for all A; > 1. Since 7^ — oo as A; — )■ oo, 
there is > 1 so that 

0<fl-^^*<ri + ::^^*<e, k>N. 



It follows that, for t > and k> N, 

I yPl"''\x, dy) < xexp | — < xe^'. (2.1.14) 

Then the desired estimate follows from (12.1.61) . (12.1.141) and Jensen's inequality. □ 

Theorem 2.1.6 Suppose that Condition \2.L2\ holds. Then for every a > we have 
Vk{t, A) — > some ft(A) uniformly on [0, a]^ as k ^ 00 and the limit function solves the 
integral equation 

v^{X) = X- f (j){vs{X))ds, X,t>0. (2.1.15) 
Jo 

Proof. The following argument is a modification of that of Aliev and Shchurenkov (1982) 
and Aliev (1985). For any n > we may write 

^oggr\^-'^') = log k(^7fc(e"'/'=))^7fc"(e-^/'=)-^] +log^7,"(e-^/'=) 
= {k^kr'Gki - khggUe-'/')) + loggUe-'^'), 

where 

Gk{z) = k^klog [^?fc(e-^/'=)e^/'=]. 
From this and (12.1.71) it follows that 

Vkit + %\ A) = Vkit, A) - %^Gkivkit, A)). 

By applying the above equation to t = 0, 1/7^, . . . , ( [7^^] — 1 ) /7fc and adding the resulting 
equations we obtain 

bkt] 

Vk{t, A) = A - 5^ Ik'GkMlk'i^ - 1), A)). 

i=l 
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Then we can write 

Vk{t,X) = X + ekit,X)- [ Gk{vk{s,X))ds, (2.1.16) 

Jo 

where 

ek{t,X) = {t-%'M)Gk{vk{%'bktlX)). 

It is not hard to see 

Gk{z) = kjk\og [1 + {k^k)-'Gk{z)e'/']. 
By Condition 12. 1.2[ for any < £ < 1 we can enlarge > 1 so that 

\Gk{z) - 4>{z)\ <e, < z <ae^'',k> N. (2.1.17) 
It then follows that 

A)| < 7^^M, 0<t,A<a, (2.1.18) 



where 



Forra > A; > let 



M = l+ sup 10(^)1. 

0<z<acBa 



Kk,n{t,X)= sup \Vn{s,X) - Vk{s,X)\. 
0<s<t 



By (12.1.161) . (12.1.171) and (12.1.181) we obtain 

Kk,n{t, A) < 2(7-iM + ea) + L ! Kk,n{s, X)ds, 0<t,X<a, 

Jo 

where L = supo<s<(icSa By Gronwall's inequality, 

Kk,nit, X) < 2{%^M + ea) exp{Lt}, < t, A < a. 

Then Vk{t,X) -> some vt{X) uniformly on [0,a]^ as A; — oo for every a > 0. From 
(12.1.161) we get (12.1.151) . □ 

Theorem 2.1.7 Suppose that cj) is given by (12.1.131) . Then for any A > there is a unique 
locally bounded positive solution t h-)- vt{X) to (12.1.151) . Moreover, the solution satisfies 
the semigroup property 

Vr+t{X) = VrOVt{X) = Vr{yt{X)), r, t, A > 0. (2.1.19) 
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Proof. By Propositions 12.1.41 and 12.1.61 there is a locally bounded positive solution to 
(12.1.151) . The proof of the uniqueness of the solution is a standard application of Gron- 
wall's inequality. The relation (12.1.191) follows from the uniqueness of the solution to 
(I2TT51) . □ 



Theorem 2.1.8 Suppose that (j) is given by (I2.1.13|) . Then there is a Feller transition 
semigroup {Qt)t>o on M+ defined by 

POO 

/ e^^yQtix, dy) = e-^^*(^\ A > 0, x > 0. (2.1.20) 
Jo 

Moreover, if 3 ^ x > 0, we have Pj^''*\xk, ■) Qt{x,-) weakly. 



Proof. By Proposition 12.1.41 and Theorems 11.1.21 and 12.1.61 there is a probability kernel 
Qt(x, dy) on IR+ defined by (12.1.201) . Moreover, we have p]^''*\xk, ■) ~^ Qt{x, ■) weakly 
if Xk X. The semigroup property of the family of kernels {Qt)t>o follows from (12.1.191) 
and (12.1.201) . For A > and x > set ex{x) = e~'^^. We denote by Di the linear span of 
{ex : A > 0}. Clearly, the operator Qt preserves Di for every t > 0. By the continuity 
of t (-7- Vt{X) it is easy to show that t i— )■ Qtex{x) is continuous for A > and x > 0. 
Then t i— )■ Qtf{x) is continuous for every f E Di and x > 0. Let Co(]R+) be the space 
of continuous functions on vanishing at infinity. By the Stone-Wei erstrass theorem, 
the set Di is uniformly dense in Co(M+); see, e.g., Hewitt and Stromberg (1965, pp.98- 
99). Then each operator Qt preserves Co(M+) and t H- Qtf{x) is continuous for every 
/ G Co(M+) and x > 0. That gives the Feller property of the semigroup {Qt)t>o- □ 

A Markov process is called a continuous-state branching process (CB-process) with 
branching mechanism if it has transition semigroup {Qt)t>G defined by (12.1.201) . It is 
simple to see that 

Qt{xi + X2,-) = Qt{xi,-) * Qt{x2,-), t,a:i,X2>0, (2.1.21) 

which is called the branching property of {Qt)t>o- The family of functions {vt)t>o is 
called the cumulant semigroup of the CB-process. By Theorem 12.1.81 the process has 
a cadlag realization. Let Vt = ^([O, oo), M+) denote the space of cadlag paths from 
[0, oo) to ]R+ furnished with the Skorokhod topology. The following theorem gives an 
interpretation of the CB-process as the approximation of the GW-process. 



Theorem 2.1.9 Suppose that Condition \2.L2\ holds. Let {x{t) : t > 0} be a CB-process 
with transition semigroup {Qt)t>o defined by (12.1.201) . If -2^(0) converges to x(0) in 
distribution, then {zk{[jkt]) : ^ > 0} converges to {x{t) : t > 0} in distribution on 
D([0,oo),M+). 
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Proof. For A > and x > set e\{x) = e Let Co(M+) be the space of continuous 
functions on M+ vanishing at infinity. By Theorem I2.1.6l it is easy to show 

Urn sup |p]^'=*'eA(x) - Qtex{x)\ =0, A > 0. 

Then the Stone-Weierstrass theorem implies 

lim sup |P^V(a:) - Qtfix)\ =0, / G Co(M+). 

By Ethier and Kurtz (1986, p.226 and pp. 233-234) we conclude that {zk{['jkt]) : t > 0} 
converges to the CB-process {x{t) : t > 0} in distribution on D{[0, oo), IR+). □ 

The convergence of rescaled Galton-Watson branching processes to diffusion pro- 
cesses was first studied by Feller (1951). Jifina (1958) introduced CB-processes in both 
discrete and continuous times. Lamperti (1967a) showed that the continuous-time pro- 
cesses are weak limits of rescaled Galton-Watson branching processes. We have followed 
Aliev and Shchurenkov (1982) and Li (2006) in some of the above calculations; see also 
Li (2011). 



2.2 Simple properties of CB-processes 

In this section we prove some basic properties of CB-processes. Most of the results pre- 
sented here can be found in Grey (1974) and Li (2000). We shall follow the treatments 
in Li (201 1). Suppose that </) is a branching mechanism defined by (12.1.131) . Then a CB- 
process has transition semigroup {Qt)t>o defined by (12.1.151) and (12.1.201) . It is easy to 
see for each x > 0, the probability measure Qt{x, ■) is infinitely divisible. Then {vt)t>o 
can be expressed canonically as 

POO 

v^(\) = ht\+ (1 -e~^")/j(dM), t>0,A>0, (2.2.1) 
Jo 

where ht > and ult{du) is a finite measure on (0,oo). From (12.1.151) we see that 
1 1-^ ft (A) is first continuous and then continuously differentiable. Moreover, we have the 
backward differential equation: 

^^Vt{X) = -(P{vt{\)), vo{\) = A. (2.2.2) 

By (12.2.21) and the semigroup property VrOVt = Vr+t for r, t > we also have the forward 
differential equation 

v,(X) = -4>(X) v^iX), vo{X) = X. (2.2.3) 
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By differentiating both sides of (12.1.151) it is easy to find 

^vt{0+) = e~'\ t > 0, (2.2.4) 
which together with (12.1.201) yields 

/■oo 

/ yQt{x, dy) = xe-^\ t > 0, x > 0. (2.2.5) 
Jo 

We say the CB-process is critical, subcritical or supercritical according as 6 = 0, > or 
< 0. 

Proposition 2.2.1 For every t > the function A ^ Vt{X) is strictly increasing on [0, oo). 

Proof. By the continuity of t i-;- fj(A), for any Aq > there is to > so that Vt{\o) > 
for < t < to- Then (12.1.201) implies Qtix, {0}) < 1 for a: > and < t < tg, and so 
A H- fj(A) is strictly increasing for < t < to- By the semigroup property of {vt)t>o we 
infer X vt{\) is strictly increasing for all t > 0. □ 



Corollary 2.2.2 The transition semigroup {Qt)t>o defined by (|2. 1.201) is a Feller semi- 
group. 



Proof By Proposition IXO for t > and A > we have Vt{\) > 0. From (12.1.201) 
we see the operator Qt maps {x ^-)■ e"^"^ : A > 0} to itself. By the Stone-Weierstrass 
theorem, the linear span of {x H- e~'^^ : A > 0} is dense in Co(IR+) in the supremum 
norm. Then Qt maps (7o(]R+) to itself. The Feller property of {Qt)t>o follows by the 
continuity of t vt(\) . □ 



Proposition 2.2.3 Suppose that A > and 0(A) 7^ 0. Then the equation (j){z) = has no 
root between A and Vt{X). Moreover, we have 

(p{z)-^dz = t, t > 0. (2.2.6) 



Proof. By (|2. 1.131) we see 0(0) = and z t-^ (f){z) is a convex function. Since 0(A) 7^ 
for some A > according to the assumption, the equation 0(z) = has at most one root 
in (0, 00). Suppose that Ao > is a root of 0(z) = 0. Then (12.2.31) implies Vt(Ao) = Aq 
for all t > 0. By Proposition 12.2. 1 1 we have Vt{\) > Xq for A > Aq and < t't(A) < Aq 
for < A < Aq. Then A > and 0(A) 7^ imply there is no root of (p{z) = between A 
andt;t(A). Fromdl^ill) wegetdlill). □ 
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Proposition 2.2.4 For any t > and \>0 let Vf{\) = {d/d\)vt{\). Then we have 

v[{X) = exp { - ^ <p'{vs{X))dsY (2.2.7) 



where 



POO 

b'{z) = b + 2cz+ m(1 -e-"")m(dM). (2.2.8) 
Jo 



Proof. Based on (12.1.151) and (12.2.21) it is elementary to see that 

|^;(A) = -<p'{v,iX))v[iX) = l^l^'t(A). 

It follows that 

|[logt.;(A)] = v[{X)-'^^v[{X) = -^'{v,iX)). 

Then we have (12.2.71) since v'q(X) = 1. □ 

Since {Qt)t>o is a Feller semigroup by Corollary I2.2.2[ the CB-process has a Hunt 
realization X = x(t), Q^). Let tq := inf{s > : x{s) = 0} denote the 

extinction time of the CB-process. 

Theorem 2.2.5 For every t > the limit Vt =tliniA^oo "^tiX) exists in (0, oo]. Moreover, 
the mapping t vt is decreasing and for any t > and x > Owe have 

QAro <t} = Q.{x(t) = 0} = exp{-xvt}. (2.2.9) 

Proof. By Proposition 12. 2. 1 1 the limit Vt =tlimA-i.oo Vt{X) exists in (0, oo] for every t >0. 
For t > r > we have 

Vt =tlini Vr{Vt^r{X)) = Vriyt-^r) < ^r- (2.2.10) 

A— )-oo 

Since zero is a trap for the CB-process, we get (12.2.91) by letting A — )■ oo in (12.1.201) . □ 

For the convenience of statement of the results in the sequel, we formulate the follow- 
ing condition on the branching mechanism: 

Condition 2.2.6 There is some constant 6 > so that 



/oo 
(j){z)~^dz < oo. 
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Theorem 2.2.7 We have vt < oo for some and hence all t > if and only if Condi- 
tion \2.2.6\ holds. 

Proof. By (12.2.101) it is simple to see that vt =tliniA->oo ft (A) < oo for alH > if and 
only if this holds for some t > 0. If Condition 12 . 2 . 6 1 holds . we can let A — oo in (|2.2.6I) 
to obtain 

roo 

/ ^{z)-^dz = t (2.2.11) 

J Vt 

and hence t;^ < oo for t > 0. For the converse, suppose that Vt < oo for some t > 0. By 
(|2.2.2I) there exists some ^ > so that (f)(9) > 0, for otherwise we would have vt{\) > A, 
yielding a contradiction. Then 0(2;) > for all z > 9 by the convexity of the branching 
mechanism. As in the above we see that (12.2.1 II) still holds, so Condition l2.2.6l is satisfied. 

□ 



Theorem 2.2.8 Let v =4,limt^oo G [0, 00]. Then for any x > Owe have 

Qx{to < 00} = exp{— x-i;}. (2.2.12) 

Moreover, we have v < 00 if and only if Condition \2.2.6\ holds, and in this case v is the 
largest root of (p^z) = 0. 

Proof. The first assertion follows immediately from Theorem 12.2.51 By Theorem 12.2.71 
we have vt < 00 for some and hence all t > if and only if Condition 12.2.61 holds . This 
is clearly equivalent to v < 00. From (12.2.1 II) it is easy to see that v is the largest root of 

0(2) = 0. □ 

Corollary 2.2.9 Suppose that Condition \2.2.6\ holds. Then for any x > we have 
Qx{to < 00} = 1 if and only ifb > 0. 

Let {Ql)t>o be the restriction to (0, 00) of the semigroup {Qt)t>o- A family of cx-finite 
measures {Ht)t>o on (0, 00) is called an entrance law for {Qt)t>o if HrQI = Hr+t for all 
r, t > 0. The special case of the canonical representation (12.2.11) with ht = for alH > 
is particularly interesting. In this case, we have 

POO 

Vt{X) = (1 - e~^")/t(dM), t > 0, A > 0. (2.2.13) 
Jo 

Theorem 2.2.10 The cumulant semigroup admits representation (12.2.131) if and only if 

POO 

0'(oo) :=6 + 2c-oo+ / um{du) = oo (2.2.14) 

Jo 

with ■ 00 = by convention. If condition (12.2.141) is satisfied, then {lt)t>o is an entrance 
law for the restricted semigroup {Ql)t>Q. 
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Proof. From (12.2.81) it is clear that the limit 0'(oo) = lim^^oo always exists in 
(— oo, oo]. By (12.2.11) we have 

POO 

v[{X) = ht+ ue~^%{du), t > 0, A > 0. (2.2.15) 
Jo 

From (12X71) and (12.2.151) it follows that 

ht = v[{oo) = exp I - ^ (f)'{vs)dsy (2.2.16) 

Then ht = for any t > implies 0'(oo) = oo. For the converse, assume that 0'(oo) = 
oo. If Condition l2.2.6l holds. by Theorem l2.2.7l for every t > we have Vt < oo, so ht = 
by (12.2.11) . If Condition 12.2.61 does not hold, then t;^ = oo for t > by Theorem 12.2.71 
Then (12.2.161) implies ht = for t > 0. That proves the first assertion of the theorem. If 
{vt)t>o admits the representation (|2.2.13l) . we can use the semigroup property of {vt)t>o 
to see 

POO POO 

/ (l-e-"")/,+t(dn) = / (l-e-"^*(^))/,(dn) 
JO Jo 

poo 

lr{dx) / (1 -e-^")Q°(a;,du) 



for r, t > and A > 0. Then {lt)t>o is an entrance law for {Q1)t>o- D 



Corollary 2.2.11 If Condition 12. 2. 6\ holds, the cumulant semigroup admits the represen- 
tation (12.2.131) and t ^ Vt = lt{0, oo) is the minimal solution of the differential equation 

^^t = -0(^J, t>0 (2.2.17) 
with singular initial condition v^^ = oo. 



Proof. Under Condition 12.2.61 for every t > we have < oo by Theorem 12.2.71 
Moreover, the condition and the convexity of z i-t- (p(z) imply </''(oo) = oo. Then we 
have the representation (12.2.131) by Theorem 12.2. 101 The semigroup property of {vt)t>o 
implies vt+s = Vt{vs) for t > and s > 0. Then t ^ vt satisfies (12.2.171) . From (12.2.1 II) 
it is easy to see vq^ = oo. Using the relation vt = limA^oo ^^t(A) it is easy to show that 
any solution t ut of (12.2.171) with mo+ = oo satisfies Ut > ik for t > 0. □ 



Corollary 2.2.12 Suppose that Condition \2.2.6\ holds. Then for any t > the function 
A H- Vt{X) is strictly increasing and concave on [0, oo), and v is the largest solution of 
the equation Vt{\) = A. Moreover, we have v =^\imt^oQ'Vt{\) for < \ < v and 
V =|Umj^oo Vt{\)for \ > v. 
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Proof. By Corollary l2.2.1 ll we have the canonical representation (12.2.131) for every t > 0. 
Since A i— j- Vt{\) is strictly increasing by Proposition 12.2. 11 the measure kidu) is non- 
trivial, so A H- Vt{X) is strictly concave. The equality v = Vt{v) follows by letting s — oo 
in vt+s = Vt{vs), where vt+s < Vs. Then v is clearly the largest solution to vt{X) = A. 
When 6 > 0, we have ?; = by Theorem 12.2.81 and Corollary 12.2.91 Furthermore, since 
0(^) > 0, from (12.2.21) we see t Vti^) is decreasing, and hence l\imt-^oo Vti^) =i 
limt_j,oo Vt = 0. If b < and < A < i;, we have A < ft(A) < Vt{v) = v for all t > 0. 
Then the limit foo(A) =t linif-t-oo ^'t(A) exists. From the relation Vt{vs{\)) = Vt^sW 
we have Wf(woo(A)) = foo(A), and hence foo(A) = v since v is the unique solution to 
ft (A) = A in (0, oo). The assertion for b < and \ > v can be proved similarly. □ 

We remark that in Theorem 12 . 2 . 1 01 one usually cannot extend {lt)t>o to ^ cr-finite en- 
trance law for the semigroup {Qt)t>o on M+. For example, let us assume Condition 12.2.61 
holds and (lt)t>o is such an extension. For any 0<r<£:<twe have 

poo pco 

ltm)> g°_,(x,{0})/,,(dx) > / e-^^-^/,(dx) 
JO Jo 

poo 

= Vr - (1 - e-™'-)/,,(du) = Vr - Vr{Vt-e). 



JO 

The right-hand side tends to infinity as r — )• 0. Then lt{dx) cannot be a cr-finite measure 
on M+. 

Example 2.2.1 Suppose that there are constants c>0, 0<a<l and b so that (f){z) = 
^^i+a _j_ Then Condition 12. 2. 6l is satisfied. Let = at and 

qi{t) = b-\l-e--''), by^O. 

By solving the equation 

^vt{X) = -cvt{XY+'' - bvtiX), vo{X) = X 

we get 

e~''*A 

^t(^) = 7 TTM ' t > 0, A > 0. (2.2.18) 

Thus Vt = c~^/°e-^*g|:;(t)-^/" for t > 0. In particular, if a = 1, then (12.2.131) holds with 

e"^* r u 
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2.3 Conditional limit theorems 

Let {Qt)t>o denote the transition semigroup of the CB-process with branching mechanism 
given by (I2.1.13I) . Let {Qf)t>Q be the restriction of {Qt)t>o to (0, oo). It is easy to 
check that Q'l{x, dy) := e''^x~^yQ1{x, dy) defines a Markov semigroup on (0, oo). Let 
q^(X) = e%t{>^) andletg;(A) = {d/dX)qt{X). Recall that ^ ^ 0'(^) is defined by ([223 ■ 
From (12.2.71) we have 

q[{X)=exp[- jy',iv,{X))ds}, (2.3.1) 
where 0o(2) = (p'{z) — b. By differentiating both sides of (12.1.201) we see 

/■oo 

/ e-^yQ'iix,dy)=exp{-xvtiX)}q[iX), A > 0. (2.3.2) 
Jo 



It follows that 



e-^yQl{x,dy)=exp[-xvt{X)- </.[,(i;,(A))ds}. (2.3.3) 
Using (12.3.31) it is easy to extend {Qf)t>o to a Feller semigroup on 

Theorem 2.3.1 Suppose that (f)'{z) oo as z oo. Let {lt)t>o be defined by (12.2.131) . 
Then for t > the probability measure Qt (0, ■) is supported by (0, oo) and (0, du) = 
ue^Ht{du). 

Proof. We first note that {vt)t>o really has the representation (12.2.131) by Theorem 12. 2. 101 



Then 



and 



/■oo 

qt{X)= / (l-e-"")e^*/i(dn), (2.3.4) 

^0 



POO 

q'^(X) = / ue-^^'e^'Htidu). (2.3.5) 
Jo 

From (I2X2I) and (12331) we see Q^(0, du) = ue^Ht{du) for u > 0. □ 

Now let X = ^t, x{t), Qx) be a Hunt realization of the CB-process with the 

augmented natural a-algebras. Let tq := mi{s > : x{s) = 0} denote the extinction 
time of X. 

Theorem 2.3.2 Suppose that b > and Condition \2.2.6\ holds. Let t > and x > 0. 
Then the distribution ofx{t) under Clx{-\f + t< tq} converges as r ^ 00 to Ql{x, ■). 
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Proof. Since zero is a trap for the CB-process, for any r > we can use the Markov 
property of {x{t) : t > 0} to see 



Q4e-"^-W|r + ^<ro] 



Um 



Qx[l{r+i<ro}] 



(2.3.6) 



Recall that tv+t = ft(fr) and t'^(O) = e By Theorem 12.2.81 and Corollary 12.2.91 we 
have liuir^oo Vr = 0. Then 

lim [e-^^W |r + t < To] = lim ^ ^ 



r— >oo 



ie'*Q,[x(t)e-^"(*)]. 



That gives the desired convergence result. □ 

It is easy to see that t i— Z{t) := e**x(t) is a positive (^t) -martingale. By the theory 
of Markov processes, for each x > there is a unique probability measure Q^. on ^) 
so that 

Ql{F) = Q,[Z{t)F] (2.3.7) 

for any ^j-measurable bounded random variable F. Moreover, under this new probability 
measure {x{t) : t > 0} is a Markov process in (0, oo) with transition semigroup {Qt)t>o- 
By a modification of the proof of Theorem 12 . 3 .21 we get the following: 

Theorem 2.3.3 Suppose that b > and Condition \2.2.6\ holds. Let x > and t > 0. 
Then for any ^t-T^^cisurahle hounded random variable F we have 

Q^[F] = lim Q,[F|r + t < tq]. (2.3.8) 

r— >oo 

By the above theorem, in the critical and subcritical cases, Q^. is intuitively the law 
of {x{t) : t > 0} conditioned on large extinction times. See Lambert (2007), Li (2000, 
20 11) and Pakes (1999) for more conditional limit theorems. 



2.4 A reconstruction from excursions 



In this section, we give a reconstruction of the sample paths of the CB-process from 
excursions. Let {Qt)t>o denote the transition semigroup of the process with branching 
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mechanism given by (I2.1.13I) . Recall that {Q1)t>o is the restriction of {Qt)t>o to (0, oo). 
Let -D(0, oo) be the space of cadlag paths t ^ Wt from (0, oo) to having zero as a trap. 
Let (^/°, s^^) be the natural cr-algebras on 1^(0, oo) generated by the coordinate process. 
For any entrance law {Ht)t>o for {Ql)t>o there is a unique (j-finite measure QH(dtf ) on 
such that Q//({0}) = and 

Qniwi^ e dxi, wt^ e dx2, . . . , wt„ e dxn) 

= Ht,{dxi)Qt^_t^{xi, dxa) ■ ■ ■ Ql^^t^_^{xn-i, dx„) (2.4.1) 

for every {ti < ■ ■ ■ < t„} C (0, oo) and {xi, . . . , x„} C (0, oo). See, e.g., Getoor and 
Glover (1987) for the proof of the existence of Qh in the setting of Borel right processes. 
Roughly speaking, the above formula means that {wt : t > 0} under Qh is a Markov 
process in (0, oo) with transition semigroup {Q1)t>o and one-dimensional distributions 

iHt)t>o- 

Theorem 2.4.1 Suppose that the condition (|2.2.14l) is satisfied. Let Q(o) be the a-finite 
measure on -D(0, oo) determined by the entrance law {lt)t>o given by (|2.2.13l) . Then for 
Q(o)-a.e. w E D{0, oo) we have Wt as t ^ 0. 

Proof. Let {Qf)t>o be the transition semigroup on M+ defined by (12.3.31) . Then we have 

Qt{x, dy) = x~^e^^yQ°{x, dy) for x,y > 0. Observe that 



D{0,oo) 



poo o 

Q(o)(dt/;) = e''ykidy) = e''—vt{0+) = l. 



Then for fixed m > we can define a probability measure Q^q-, (dw) := e WuQ,{o) (dw) on 
D{0, oo). Under this measure, the coordinate process {wt : < t < m} is an immigration 
process with transition semigroup {Qt)t>o and one-dimensional distributions 

g,^(0,d2/) = eVt(d|/), 0<t<u. 

By the uniqueness of the transition law of the immigration process we have ivt — t- as 
t — )■ for Q^Q^-a.e. w G -D(0, oo). Note that Q"g^(dw) and Q(o)(dw) are absolutely 
continuous with respect to each other on Du{0, oo) := {w G D{0, oo) : Wu > 0}. Since 
D{0, oo) = {0} U (U^=iL'i/„(0, oo)) and Q(o)({0}) = 0, we have ^ as t ^ for 
Q(o)-a.e. If G -D(0, oo). □ 

By Theorem l2.4.1l we can regard Q(o) as a a-finite measure on the space of excursions 
Do[0,oo) := {w G D{0,oo) : wq = 0}. We call Q(o) an excursion law for the CB- 
process. Using the excursions, we can give a reconstruction of the sample paths of the 
CB-process. Let x > and let N(dw) be a Poisson random measure on Dq[0, oo) with 
intensity xQ(o)(dw). We define the process {Xf : t > 0} by Xq = x and 



Xt= / WtN{dw), t>0. (2.4.2) 

'A)[0,oo) 
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Theorem 2.4.2 For t > let be the a -algebra generated by the collection of random 
variables {N{A) : A e ^°}. Then {(X^, : t > is a realization of the CB-process. 

Proof We first remark that the random variable Xt has distribution ■) on IR^.. In 
fact, for any t > and A > we have 

P[exp{-AXi}] = exp { - X / (1 - e-^'^OQwldu;)) 

JDo[0,oo) 

= expj -x^ (1 -e-^")/t(dz)} = e^^{-xvt{\)}. 

Let t > r > and let /i be a bounded positive function on -Do[0, oo) measurable relative 
to . For any A > we have 



exp 



h{w)N{dw) - \Xt 



A)[0,oo) 



-h{w)—Xwt 



)Q(0)(d 



exp \ — X (l - 

Jdo[0,oo) 

exp{-x / (l-e-'^("'))Q(o)(di/;)| 

^DnfO.oo) 



• exp ^ — X j e 



-M-)(i_e-^-*)Q(o)(dw: 



where we made the convention e °° = 0. By the Markov property of Q(o) we have 
I e-'^("')(l-e-^"'*)Q(o)(dM;) 

JDo[0,oo) 



e-'^("')Q(o)(dw;) / (l-e-^^)Q°_,K,dz) 



DoIO.oo) 



(1 





-tOrt't-r{A) 



)Q(o)(d 



w). 



It follows that 



exp 



h{w)N{dw) - \Xt 



Do[0,cx>) 



exp 



{-/ 



-'^('^)o-'^'-'^t-r-(-^) 



)Q(o)(d 



w] 



exp 



/i(tf;)A^(dtf;) -X,t;t„^(A) 



'Do[0,oo) 

Then {{Xt, 5^t) : t > 0} is a Markov process with transition semigroup {Qt)t>o- 



□ 
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Chapter 3 

Structures of independent immigration 



In this chapter we study independent immigration structures associated with CB-processes. 
We first give a formulation of the structures using skew convolution semigroups. Those 
semigroups are in one-to-one correspondence with infinitely divisible distributions on M+. 
We show the corresponding immigration process arise as scaling limits of Galton-Watson 
processes with immigration. We discuss briefly limit theorems and stationary distribu- 
tions of the immigration superprocesses. The trajectories of the immigration processes 
are constructed using stochastic integrals with respect to Poisson random measures deter- 
mined by entrance laws. 

3.1 Formulation of immigration processes 

In this section, we introduce a generalization of the CB-process. Let {Qt)t>o be the tran- 
sition semigroup defined by (12.1.151) and (12.1.201) . Let {'jt)t>o be a family of probability 
measures on M+. We call (7t)t>o a skew convolution semigroup (SC-semigroup) associ- 
ated with {Qt)t>o provided 



-fr+t = {irQt) * It 



r,t > 0. 



(3.1.1) 



It is easy to show that (13.1.11) holds if and only if 



Ur+t{,X) = Ut{X) + Ur{Vt{X)), 



r,t, A > 



(3.1.2) 



where 




(3.1.3) 



The concept of SC-semigroup is of interest because of the following: 
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Theorem 3.1.1 The family of probability measures {jt)t>o on R+ is an SC-semigroup if 
and only if 

Q]{x,-) ■.= Qt{x,-)*7t, t,x>0 (3.1.4) 
defines a Markov semigroup {Q])t>o on M+. 

Proof. Let (7t)t>o probability measures on IR+ and let Q]{x, ■) be the probability kernel 
defined by (13.1.41) . Then we have 



oo 



e-y^QUx,dy) = exp{-xvt{X)-Ut{X)}, t,x,A>0. (3.1.5) 

Using this relation it is easy to show that {Q])t>o satisfies the Chapman-Kolmogorov 
equation if and only if (13.1.21) is satisfied. That proves the result. □ 

If {y{t) : t > 0} is a positive Markov process with transition semigroup {Qj)t>o given 
by (13.1.41) . we call it an immigration process or a CBI-process associated with {Qt)t>Q- 
The intuitive meaning of the model is clear in view of (13.1.11) and (13.1.41) . From (13.1.41) 
we see that the population at any time t > is made up of two parts; the native part 
generated by the mass x > has distribution Qt{x,-) and the immigration in the time 
interval (0, t] gives the distribution 7^. In a similar way, the equation (13.1.11) decomposes 
the mass immigrating to the population during the time interval (0, r + t] into two parts; 
the immigration in the interval (r, r + t] gives the distribution 7^ while the immigration 
in the interval (0, r] generates the distribution 7^ at time r and gives the distribution 'yrQt 
at time r + t. It is not hard to understand that (13.1.41) gives a general formulation of the 
immigration independent of the state of the population. 

Theorem 3.1.2 The family of probability measures {'yt)t>o on M+ is an SC-semigroup if 
and only if there exists ip E so that 

^ e-^y^tidy) = exp I - ^ ^(t;,(A))ds}, t, A > 0. (3.1.6) 

Proof. It is easy to check that for any ^ e ^ the family (7t)t>o defined by (13.1.61) is 
an SC-semigroup. Conversely, suppose that (7t)j>o is an SC-semigroup. For t, A > let 
Ut{X) be defined by (13.1.31) . Then t Ut{X) is increasing. By Lebesgue's theorem, the 
limit 

i^tW := lim s~^[mj+s(A) - Mt(A)] = lim s^'^UsivtiX)) 

exists for almost all t > 0; see, e.g., Hewitt and Stromberg (1965, p. 264). By the con- 
tinuity of t I— 7- Vt{X), there is a dense subset D of (0, 00) so that the following limits 
exist: 

ijoW ■= lim s~\s{X) = lim - e'^^^^^], XeD. (3.1.7) 
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For u G [0, oo] and A G (0, oo) let ^{u, A) be defined by (I1.2.4I) . Then ^{u, A) is jointly 
continuous in (m, A). By (13.1.71) we have 

V'o(A) = lim s-\s{\) = lim s"^ / ^{u,\)Gs{du) (3.1.8) 

«^0+ «^0+ i[0,oo] 

for X E D, where Gs{du) = (1 — e~")7s(du). Then for some 5 > we have 

sup s~^G<i([0, oo]) < oo, 

0<s<<5 

so the family of finite measures {s^^Gg : < s < 5} on [0, oo] is relatively compact. 
Suppose that s„ — > and s~^Gs„ some G weakly as ?i — oo. Then 

lim s-^UsM = lim f ^{u, A)a„(dM) = [ ^u, X)G{du), A > 0. 

""^"^ J[0,oo] J[0,oo] 

Thus we can extend to a continuous function on (0, oo) given by 



^o(A) = / ^{u,\)Gidu) = Gi{oo}) + h\+ (1 -e-"^)/(rfM), 

i[0,oo] i(0,oo) 

where h = G{{0}) and l{du) = (1 — e~")~^Gs(dn). By a standard argument one sees 
that (13X81) holds actually for all A > 0. From (l3T2l) and (l3X8l) it follows that 

I^+«s(A)|^^, = D+w,K(A))|^^^ = ^oK(A)), 

where D"^ denotes the right derivative relative to s > 0. Here the right-hand side is 
continuous in t > 0. Thus 1 1— )■ Ut{X) is continuously differentiable and (13.1.61) holds with 
^Ij(X) = ^^o(A) for A > 0. By letting A 0+ in (l3T6l) one sees G({oo}) = 0. Then 

ij = iJoey. □ 

By Theorem 13.1.21 there is a 1-1 correspondence between SC-semigroups and in- 
finitely divisible distributions on M+. Then the theorem generalizes the 1-1 correspon- 
dence between classical convolution semigroups and infinitely divisible distributions. In 
fact, from (13.1.11) it is easy to see that (7j)t>o reduces to a classical convolution semigroup 
if Qt is the identity operator for all t > 0. As a consequence of Theorems ll.2.4l and l3.1.2[ 
an SC-semigroup {'jt)t>o always consists of infinitely divisible distributions. 

Now let us consider a transition semigroup {Qj)t>o defined by (13.1.41) with the SC- 
semigroup (7t)t>o given by (13.1.61) . If an immigration process has transition semigroup 
{Qt )t>o, we say it has branching mechanism (p and immigration mechanism ip. It is easy 
to see that 

^ e-'yQ]{x,dy) = exp\-xvt{X) - f ij{vs{X))ds}. (3.1.9) 
Jo ■' 
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Recall that the function ip ^ has the representation 

poo 

^p{z)=l3z+l (1 -e~"")n(du), ^ > 0, (3.1.10) 
Jo 

where /3 > is a constant and (1 A u)n{(iu) is a finite measure on (0, oo). By (13.1.91) and 
(12.2.41) one can show 

/•oo rt 

/ yQj{x,dy) = xe-'' -^'{^) e-'^ds, (3.1.11) 
Jo Jo 

where 

/•oo 

^'(0) = 13+ un{du). (3.1.12) 
Jo 

The following theorem gives a necessary and sufficient condition for the ergodicity of 
the semigroup {Q])t>o- 

Theorem 3.1.3 Suppose that b > and (f){z) ^ for every z > 0. Then Q]{x,-) 
converges to a probability measure r] on M+ as t ^ oo if and only if 

^ il)(z) 

-^dz < oo for some X > 0. (3.1.13) 

(f){z) 

(13. 1.131) holds, the Laplace transform off] is given by 

L^(A)=exp|-^ i){vs{\))dsy A > 0. (3.1.14) 

Proof. Since (j){z) > for all ^ > 0, from (12.2.21) we see t i-)- Vt{X) is decreasing. Then 
(|2.2.6I) implies limj_^oo vt{X) = 0. By (13.1.91) we have 

Um^ e~^yQ]{x,dy) = exp\^- tlj{vs{X))ds'j (3.1.15) 



for every A > 0. A further application of (|2.2.2I) gives 

r-i /•A 

■tjj' 



2lj{Vs{X))ds = I '^dz. 



It follows that 



MX) 



^{v,{X))ds= l^^dz, 



rv^) 

which is a continuous function of A > if and only if (13.1.131) holds. Then the result 
follows by (13.1.151) and Theorem [TX2l □ 
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Corollary 3.1.4 Suppose that b > 0. Then Q]{x, ■) converges to a probability measure rj 
on ]R+ as t ^ oc if and only if \ogun{du) < oo. In this case, the Laplace transform 
off] is given by (I3.1.14I) . 

Proof. We have (p^z) = bz + o{z) as z — )■ 0. Thus (13.1.131) holds if and only if 



-dz < oo for some A > 0, 







z 



which is equivalent to 



^ dz f°° f^^ 1 — e~y 
— / (l - e"'''")n(dM) = / n{du) I dy < oo 







for some A > 0. The latter holds if and only if logMn(dM) < oo. Then we have the 
result by Theorem [3TT3J □ 

In the situation of Theorem 13. 1.31 it is easy to show that 77 is a stationary distribution 
for {Ql)t>Q- The fact that the CBI-process may have a non-trivial stationary distribution 
makes it a more interesting model in many respects than the CB-process without immi- 
gration. Note also that the transition semigroup {Q\)t>o given by (12.3.31) is a special case 
of the one defined by ( 13.1.91 ). 

Theorem 3.1.5 Suppose that b > and let g{(A) be defined by (12.3.11) . Then for every 
\> the limit q'{X) ■.=l\imt^oo QtW exists and is given by 

g'(A)=exp|-^ <j)',{vs{X))dsY A > 0. (3.1.16) 

Moreover, we have ^'(0+) = g'(0) = 1 if and only if ulogum^du) < 00. The last 
condition is also equivalent to q'{X) > Ofor some and hence all A > 0. 

Proof. The first assertion is easy in view of (12.3.11) . By Corollary 13.1.41 we have 
u\ogum(du) < 00 if and only if A 1— ^'(A) is the Laplace transform of a proba- 
bility r] on M+. Then the other two assertions hold obviously. □ 

Theorem 3.1.6 Suppose that b > and (p'{z) 00 as z ^ 00. Then Q\{x, ■) converges 
as t — 7- 00 to a probability r] on (0, 00) if and only if u\ogum{du) < 00. If the 
condition holds, then r] has Laplace transform = q' given by (I3.1.16|) . 

Proof. By Corollary 13.1.41 and Theorem 13.1.51 we have the results with r] being a prob- 
ability measure on IR+. By Theorem 12.3. II the measure Qt{0, ■) is supported by (0, 00), 
hence Qt{x, ■) is supported by (0, 00) for every x > 0. From (12.3.51) we have L^(oo) < 
g;(oo) = for t > 0. That implies r/({0}) = 0. □ 
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Example 3.1.1 Suppose that c>0, 0<a<l and b are constants and let 0(2) = 
_|_ for 2; > 0. In this case the cumulant semigroup {vt)t>o is given by (I2.2.18I) . 
Let l3 > and let ip{z) = /3z°' for z > 0. We can use (13.1.91) to define the transition 
semigroup {Q])t>o- It is easy to show that 



The concept of SC-semigroup associated with branching processes was introduced in 
Li (1995/6, 1996). Theorem 13. 1 .2l can be regarded as a special form of main theorem of Li 
(1995/6). Theorem 13. 1 .31 and Corollary 13 . 1 .41 were given in Pinsky (1972). Other results 
in this section can be found in Li (2000). 

3.2 Stationary immigration distributions 

In this section, we give a brief discussion of the structures of stationary distributions of 
the CBI- processes. The results here were first given in Li (2002) in the setting of measure- 
valued processes. Given two probability measures 771 and r/2 on M+, we write ?7i ^ 772 if 
r/i * 7 = ?72 for another probability measure 7 on M+. Clearly, the measure 7 is unique if 
it exists. Let {Qt)t>o be the transition semigroup defined by (|2. 1.151) and (12.1.201) . where 
{vt)t>o has the representation (|2.2.1I) . Let {Q°)t>o be the restriction of {Qt)t>o to (0, 00). 
Let (S'*{Q) denote the set of probabilities t] on M+ satisfying rjQt ^ r] for all t > 0. 

Theorem 3.2.1 For each r] G (S'*{Q) there is a unique SC-semigroup {pit)t>o associated 
with {Qt)t>o such that rjQt * = vfaf t > and 1] = limt^oo It- 
Proof. By the definition of (S'*(Q), for each t > there is a unique probability measure 
7t on M+ satisfying 77 = {'qQt) * 7*. By the branching property of {Qt)t>o one can show 
(yUi * ^2)Qt = il^iQt) * (yU2Qi) for any t > and any probability measures and ^2 on 
]R+. Then for r, t > we have 

{riQr+t) * Ir+t = {vQt) *lt = {[{'nQr) * lr]Qt} * It = (vQr+t) * {irQt) * It- 

A cancelation gives (13.1.11) . so (7j)j>o is an SC-semigroup associated with {Qt)t>o- Now 
for every A > the function t 1— )■ L^t{X) is decreasing. By the relation 77 = (i^Qt) * 7t 
one can see t t-^ LriQ^{X) is increasing. Then there are probability measures ?7j and ?7p on 
M_i_ so that rii*rip = 7] and 




A > 0. 



(3.1.17) 



L^^(A) = lim L^Q,(A), L 



(A) = lim L,,(A). 
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These imply rji = limt^oo vQt and rjp = limj^oo7t- It follows that rji is a stationary 
distribution of {Qt)t>o, so we must have rji = 60 and rjp = rj. □ 

In the situation of Theorem l3.2.1l it is easy to see the measure t] E S'*(Q) is the unique 
stationary distribution of the transition semigroup {Qt )t>o defined by (13. 1.41) . Then we 
can identify (S'*{Q) with the set of stationary distributions of immigration processes asso- 
ciated with {Qt)t>Q- As a consequence of Theorem l3.1.3l and l3.2.ll every /i G <S'*{Q) is 
infinitely divisible. Recall that we write /i = I{h, I) if fi is an infinitely divisible proba- 
bility measure on M+ with ^ := — logL^ given by (11.2.31) . Let S{Q°) denote the set of 
excessive measures u for {Ql)t>o satisfying 

/ (1 A u)u{&a) < 00. 
Jq 

The following result gives some characterizations of the set (S'*{Q). 

Theorem 3.2.2 Let rj = I{I3, v) be an infinitely divisible probability measure on M+. 
Then rj G (S'*{Q) if and only if{f3, u) satisfy 

(3ht < (3 and (3lt + uQl < u, t>0. (3.2.1) 

In particular, ifu G S{Q°), then 77 = /(O, v) G S*{Q). 

Proof. It is easy to show that rjQt = I{(3t,i't), where (3t = (3ht and ut = (3lt + yQ^. Then 
TlQt ^ V holds if and only if (13.2.11) is satisfied. The second assertion is immediate. □ 



3.3 Scaling limits of discrete immigration models 

In this section, we prove a limit theorem of rescaled Galton- Watson branching processes 
with immigration, which leads to the CBI-processes. This kind of limit theorems were 
studied in Aliev and Shchurenkov (1982), Kawazu and Watanabe (1971) and Li (2006) 
among many others. 

Let g and h be two probability generating functions. Suppose that {^n,i '■ n,i = 
1, 2, . . .} and {?7„ : n = 1, 2, . . .} are independent families of positive integer-valued 
i.i.d. random variables with distributions given by g and h, respectively. Given another 
positive integer-valued random variable y{0) independent of {Cn,i} and {i]n}, we define 
inductively 



j/(n-l) 
i=l 



(3.3.1) 
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Then {y{n) : n = 0, 1, 2, . . .} is a discrete-time positive integer- valued Markov chain 
with transition matrix Q{i,j) determined by 

oo 

^Q(^,J>'■ = 9{zyh{z), \z\ < 1. (3.3.2) 

The random variable y{n) can be thought of as the number of individuals in generation 
n > of an evolving particle system. After one unit time, each of the y{n) particles splits 
independently of others into a random number of offspring according to the distribution 
given by g and a random number of immigrants are added to the system according to 
the probability law given by h. The n-step transition matrix Q"'{i,j) of {y{n) : n = 
0, 1,2,.. .} is given by 

oo n 
3=0 j=l 

where g'^{z) is defined by g'^{z) = g{g"'~^{z)) successively with g^{z) = z. We call 
any positive integer- valued Markov chain with transition probabilities given by (13.3.21) or 
(13.3.31) a Galton-Watson branching process with immigration (GWI-process) with param- 
eters {g,h). When h = 1, this reduces to the GW-process defined in the first section. 

Suppose that for each integer A; > 1 we have a GWI-process {yk{n) : n > 0} with 
parameters {gk,hk). Let Zk{n) = yk{n)/k. Then {zk{n) : > 0} is a Markov chain 
with state space := {0, 1/k, 2/k, . . .} and ra-step transition probability Q^{x, dy) de- 
termined by 

« n 

/ e-'yQUx,dy)=gUe^'/'')''^l[h{gi-\e^'/')), A > 0. (3.3.4) 

Suppose that {'jk} is a positive real sequence so that 7^ — ?■ 00 increasingly as A; — )■ 00. 
Let [jkt] denote the integer part of 7fct > 0. In view of (13.3.41) . given 2:^(0) = x the 
conditional distribution ■) of 2fc([7fct]) on Ek is determined by 

/ e"'yQi'\x,dy) 

J Ek 

= exp I - xffc(t, A) - y Hk{vk{s,\))dsY (3.3.5) 

where Vk{t,X) is given by (12.1.71) and 

Hk{X) = -lklogh{e-'/''), A>0. 
For any 2; > let Gk{z) be defined by (12.1.81) and let 

Hk{z) = - hk{e-'/% (3.3.6) 
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Condition 3.3.1 There is a function ip on [0, oo) such that Hk{z) ip{z) uniformly on 
[0, a] for every a > as k ^ oc. 

It is simple to see that Hk E . By Theorem ll.2.2[ if the above condition is satisfied, 
the limit function ^ has the representation (I3.1.10I) . A different of proof of the following 
theorem was given in Li (2006). 

Theorem 3.3.2 Suppose that Conditions \2.L2\ and \3.3.1\ are satisfied. Let {y{t) : t > 0} 
be a CBI-process with transition semigroup {Qt)t>o defined by (I3.1.9I ). If z^if)) converges 
to y{0) in distribution, then {^^([7^^]) : t >0} converges to {y(t) : t > 0} in distribution 
on L)([0,oo),M+). 

Proof. By Theorem 12. 1 .61 for every a > we have Vk{t, A) — )■ ft(A) uniformly on [0, a]^ 
as — J- 00. For A > and a; > set ex{x) = e~^^. In view of (13.3.51) we have 

lim sup \Q^^''^^ex{x) - Q]ex{x) \ = 
for every t > 0. Then the result follows as in the proof of Theorem 12. 1.91 □ 

Example 3.3.1 In a special case, we can give a characterization for the CBI-process in 
terms of a stochastic differential equation. Let m = E[^i 1]. From (13.3.11) we have 

y{n) - y{n - 1) = ^/y{n - 1) V - (1 - m)y{n - 1) + r]n. 

Then it is natural to expect that a typical CBI-process would solve the stochastic differen- 
tia equation 

dy{t) = ^2cy{t)dB{t) - by{t)dt + /3dt, t > 0, (3.3.7) 

where {B{t) : t > 0} is a Brownian motion. The above equation has a unique positive 
strong solution; see Ikeda and Watanabe (1989, pp. 235-236). In fact, the solution {y{t) : 
t > 0} has transition semigroup given by (13.1.171) with a = 1. Let C^(M+) denote the set 
of bounded continuous real functions on IR+ with bounded continuous derivatives up to 
the second order. Then {y{t) : t >0} has generator A given by 

In particular, for /3 = the solution of (13.3.71) is called Feller's branching diffusion. 
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3.4 A reconstruction of the sample paths 

In this section, we give a reconstruction of the sample paths of the CBI-process using a 
Poisson random measure. Let us consider an entrance law {Ht)t>o for {Qt)t>o and let 
be the cr-finite measure on D{0, oo) determined by (12.4.11) . Suppose that N{ds, dw) is a 
Poisson random measure on (0, oo) x D(0, oo) with intensity dsQ//(dtf ). We define the 
process 

Yt= [ [ Wt-sN{ds,dw), t>0. (3.4.1) 

J{0,t) Jd{0,od) 

A special form of the following theorem was established in Pitman and Yor (1982): 



Theorem 3.4.1 For t > let be the a-algebra generated by the collection of random 
variables {N{{Q,u] x A) : A e £^t\,0 < u < t}. Then {ijt,'^^ : ^ > 0} is an 
immigration process with transition semigroup {Q^)t>o given by 

e-^yQf{x,dy)=exp[-xvt{X) - j^s j^^ {1 - e-'y)H,{dy)y (3.4.2) 



Proof. It is easy to show that Yt has distribution Q^{0, ■) on M_|_. Let t > r > n > and 
let hhe a bounded positive function on D(0, oo) measurable relative to £^r~u- A > 
we can see as in the proof of Theorem 12.4.21 that 



exp 



^D{0,oo) 
t 



^D{0,oo) 



h{w)N{ds,dw) - XYt 

[h{w)l{s<u} + Xwt-s]N{ds,dw) 
expj-/" ds [ (1 -e-'^("')^{-^">e-^""-^)Qj^(dM;)| 

Jo J D {0,00) ^ 

expj-/" ds I {I - e-^^'"h-'''~^^^^'"^-^)ClH{dw) 
[ ds [ (1 -e"'''-'-^^)"'— )QH(di/;)| 

Ju iD(0,oo) 



■ exp 

■ exp 
exp 

• exp 



ds 



(1 



-\wt- 



D{0,oo) 

h{w)N{ds, dw) - YrVt-r{X) 



JD{0,oo) 



j\s j^^{l-e~'y)H,^,{dy)y 
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where we have used the Markov property (12.4.11) for the third equality. That shows 
{{Yt, '^t) : t > 0} is a Markov process in M+ with transition semigroup {Q^)t>o- O 

Let Px{dw) denote the distribution on D[0, oo) of the CB-process {x{t) : t > 0} with 
x{Q) = X. Suppose that ip e y he given by (13.1.101) . If the condition (12.2.141) is satisfied, 
we can define a a-finite measure Qnidw) on D[0, oo) by 

/■oo 

Qnidw) = /3Q(o)(dii;) + / n(dx)P,(dw;). (3.4.3) 

This corresponds to the entrance law {Ht)t>Q for {Qf)t>o defined by 

Ht = (3k + / r2(dx)Qt(x,-), t>0. (3.4.4) 
Jo 

In this case, it is easy to show that 

/■oo 

/ {l-e-^y)Ht{dy) = ^{vt{X)), t > 0, A > 0. (3.4.5) 
Jo 

Then from Theorem 13 .4. II we obtain 

Corollary 3.4.2 Suppose that (12.2.141) is satisfied and {Ht)t>o is given by (13.4.31) . Then 
{{Yt, ^t) ■ t > 0} is an immigration process with transition semigroup {Q])t>o given by 
(I3T91) . 
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Chapter 4 

Martingale problems and stochastic 
equations 



Martingale problems play a very important role in the study of Markov processes. In this 
chapter we prove the equivalence of a number of martingale problems for CBI-processes. 
From the martingale problems we derive some stochastic equations. Using the stochastic 
equations, we give simple proofs of Lamperti's transformations on CB-processes and 
spectrally positive Levy processes. 



4.1 Martingale problem formulations 

In this section we give some formulations of the CBI-process in terms of martingale prob- 
lems and prove their equivalence. Suppose that (0, tp) are given respectively by (12.1.131) 
and (|3. 1.101) with un{du) being a finite measure on (0, oo). For / G C^(M+) define 



Lf{x) = cxf'ix) + x [f{x + z)- f{x) - zf{x)]m{dz) 

Jo 

poo 

+ {/3-bx)f{x)+ / [f{x + z)-f{x)]n{dz). (4.1.1) 



We shall identify the operator L as the generator of the CBI-process. For this purpose we 
need the following: 

Proposition 4.1.1 Let {Q2)t>o be the transition semigroup defined by (12.1.201) and (13.1.91) . 

Then for any t > and A > we have 

/ e-^^Q7(x, dy) = e-^^ + ds [x0(A) - V^(A)]e-^^Q](x, dy). (4.1.2) 
Jo Jo Jo 
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Proof. Recall that ^{(A) = {d/d\)vt{\). By differentiating both sides of (13.1.91) we get 



ye-y'QUx,dy) 



From this and (|2.2.3I) it follows that 



d_ 

di 



d 



e-y^Q2{x, dy) = - [x^^vtiX) + ^(t;*(A))J j e'^^Q^x, dy) 

xm<W-i^W / e-y'QUx,dy) 
^ Jo 

ij'{vs{\))-^^vs{\)ds / e-y'Q:{x,dy) 



xmv'tW-i'W / e-y'Q:{x,dy) 
^ Jo 

POO 

+ 0(A) / i;'ivs{\)Ki\)ds / e~y'Q:ix,dy) 
Jo Jo 

[x<P{\)-i:{X)]e~y'Q:{x,dy). 



That gives (l4T2l) . 



□ 



Suppose that (Q, ^, ^^j, P) is a filtered probability space satisfying the usual hypothe- 
ses and {y{t) : t > 0} is a cadlag process in R_|_ that is adapted to i^t)t>o and satisfies 
P[?/(0)] < oo. Let us consider the following properties: 



(1) For every T > and A > 0, 

exp I - VT--tWy{t) 



T-t 



0<t<T, 



is a martingale. 
(2) For every A > 0, 



H,iX) := exp { - Xy{t) + /"^[^(A) - y{s)<f){X)]d. 



t > 0, 



is a local martingale. 

(3) (a) The process {y{t) : t >0} has no negative jumps. Let A^(ds, dz) be the optional 
random measure on (0, oo)^ defined by 

N{ds,dz) = ^ l{Ay{s)=io}hsAy{s))ids, dz), 

s>0 
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where Ay{s) = y{s) — y{s—), and let iV(ds, dz) denote the predictable compen- 
sator of N(ds, dz). Then iV(ds, d^;) = y{s—)dsm(dz) + dsn{dz). 

(b) If we let N{ds, dz) = N{ds, dz) - N{ds, dz), then 



ft 

rd 



y{t) = 2/(0) + M^ + M^+ 13+ zn{dz) - hy{s) 







ds, 



where t n- M^ is a continuous local martingale with quadratic variation 2cy{t—)dt 
= 2cy{t)dt and 



pt poo 

t^Mf= / zN{ds,dz) 
Jo Jo 



is a purely discontinuous local martingale. 
(4) For every / G C^(M+) we have 

f{y{t)) = f{y{0)) + /* Lf{y{s))ds + local mart. 

Theorem 4.1.2 The above properties (1), (2), (3) and (4) are equivalent to each other. 
Those properties hold if and only if {{y(t),^t) t > 0} is a CBI-process with parameters 

Proof. Clearly, (1) holds if and only if {y{t) : t > 0} is a Markov process relative to 
(^i)i>o with transition semigroup {Qt)t>o defined by (13.1.91) . Then we only need to prove 
the equivalence of the four properties. 

(1)^(2): Suppose that (1) holds. Then {y{t) : t > 0} is a CBI-process with transition 
semigroup {Q])t>o given by (13.1.91 ). By (14.1.21) and the Markov property it is easy to see 
that 



Yt{X) := e-^^W + / [^(^) _ y{s)(j){X)]e~^y^'^ds 
Jo 

is a martingale. By integration by parts applied to 

Zt{X) := e-^^W and Wt{X) := exp | ^ [^(A) - y{s)<P{X)]ds'j (4.1.3) 

we obtain 

dHtiX) = e~^y^'-^dWtiX) + WtiX)de-^y^'^ = WtiX)dYtiX). 
Then {Ht{X)} is a local martingale. 
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(2) ^(3): For any A > define Zt{X) and Wt{X) by (14X31) . We have Zt{X) = 
Ht{X)Wti\)~^ and so 

dZt(A) = Wt{X)-'dHt{X) - Zt4X)[^{X) - y{t-)<f){X)]dt (4.1.4) 

by integration by parts. Then {Zt{X)} is a special semi-martingale; see, e.g., Dellacherie 
and Meyer (1982, p. 213). By Ito's formula we find the {y{t)} is also a special semi- 
martingale. We define the optional random measure N{ds, dz) on [0, oo) x M by 

N{ds,dz) = l{Ay(s)^o}S(s,Ay(s)){ds, dz), 

where A?/(s) = y{s) — y{s—). Let N{ds,dz) denote the predictable compensator of 
N{ds, dz) and let N{ds, dz) denote the compensated random measure; see Dellacherie 
and Meyer (1982, pp.371-374). It follows that 

y{t) = y{0) + Ut + M,' + Mt, (4.1.5) 

where {Ut} is a predictable process with locally bounded variations, {Mf } is a continuous 
local martingale and 

Mf=[ [ zN{ds,dz), t>0, (4.1.6) 
Jo Jr 

is a purely discontinuous local martingale; see Dellacherie and Meyer (1982, p. 353 and 
p. 376) or Jacod and Shiryaev (2003, p. 84). Let {Ct} denote the quadratic variation pro- 
cess of {M^}. By Ito's formula, 

Zt{X) = Zo{X) - X f Z,4X)dU, + ^A^ fz,_{X)dCs 
^ Jo ^ Jo 

+ / / (A) (e^^^ - 1 + 2A)iV(ds,d^) + local mart. (4.1.7) 
Jo Jr 

In view of (14.1.41) and (14.1.71) we get 

[l/(t)0(A) - 'ipiX)]dt = Ix'^dCt - XdUt + / (e""^ - 1 + zX)N{dt, dz) 

2 Jr 

by the uniqueness of canonical decompositions of special semi-martingales; see Del- 
lacherie and Meyer (1982, p. 213). By substituting the representation (12.1.131) of cf) into 
the above equation and comparing both sides it is easy to find that (3. a) and (3.b) hold. 

(3) ^(4): This follows by Ito's formula. 

(4) ^(1): Let G = G{t, x) G ^^'^([O, T] x M+). For < t < T and A; > 1 we have 

oo 

G{t, y{t)) = G{0, y{0)) + ^ [G{t A j/k, y{t A (j + l)/k)) - G{t A j/k, y{t A j/k))] 

j=0 
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+ J] [Git A (j + l)/k, y{t A (j + l)/k)) - G{t A j'/fc, yit A (j + 1)/A:))] , 
j=0 

where the summations only consist of finitely many non-trivial terms. By applying (4) 
term by term we obtain 



G{t,y{t))=G{0,y{0)) + Y, 

j=0 



l'tAij+l)/k 



thj/k 



[P-hy{s)]G'{tA3lKy{s)) 



+ cy{s)G'Ut A 3/k, y{s)) + y{s) / G{t A j/k, y{s) + z) 

Jo 

- Git A 3/k, yis)) - zG'yit A 3/k, m(dz) 

[^(t A j-A, + ^) - G(t A yis))] nidz)]ds 



i=o 



tA{j + l)/fe 



G;(s,2/(tA(j + l)A))ds + Mfc(t), 



where {Mfe(t)} is a local martingale. Since {yit)} is a cadlag process, letting A; — oo in 
the equation above gives 



= G(0,y(0)) + j \G',is,yis))-hyis)G'yisMs)) 

^ /"OO 

+ cyis)Gl^is, yis)) + yis) / + z) 

Jo ^ 

~Gis,yis)) - zG'^is,yis)) m(dz) 

Gis,yis) + z)- Gis,yis)) n(d^)|ds + M(t), 

where {M(t)} is a local martingale. For any T > and A > we may apply the above to 



Git, x) = exp I — VT~t{X)x — J 



i^{vsW)ds 



to see t Git, yit)) is a local martingale. 



□ 



The above property (4) implies that the generator of the CBI-process is the closure of 
the generator L in the sense of Ethier and Kurtz (1986). This explicit form of the generator 
was first given in Kawazu and Watanabe (1971). 



4.2 Stochastic equations of CBI-processes 

In this section we establish some stochastic equations for the CBI-processes. The reader 
may refer to Dawson and Li (2006, 2010) and Fu and Li (2010) for more results on this 



48 CHAPTER 4. MARTINGALE PROBLEMS AND STOCHASTIC EQUATIONS 



topic. Suppose that {(f), tp) are given respectively by (12.1.131) and (13.1.101) with un{du) 
being a finite measure on (0, oo). Let {Qt)t>o be the transition semigroup defined by 
(12.1.201) and (13.1.91 ). In this section, we derive some stochastic equations for the CBI- 
processes. 

Suppose that (fi, P) is a filtered probability space satisfying the usual hy- 

potheses. Let {B{t) : t > 0} be an (^j)-Brownian motion and let {po{t) : t > 0} and 
{pi{t) : t > 0} be (^t)-Poisson point processes on (0, oo)^ with characteristic measures 
m{dz)du and n(dz)du, respectively. We assume that the white noise and the Poisson pro- 
cesses are independent of each other. Let No^ds, dz, du) and Ni(ds, dz, du) denote the 
Poisson random measures on (0, oo)^ associated with {poit)} and {pi{t)}, respectively. 
Let No{ds, dz, du) denote the compensated measure of N^^ds, dz, du). Let us consider 
the stochastic integral equation 



where A''o(ds, dz, du) = A''o(ds, dz, du) — ds'm{dz)du. We understand the last term on 
the right-hand side as an integral over the set {(s, 2;, m) : < s < t, < 2 < oo, < m < 
y(s—)} and give similar interpretations for other integrals with respect to Poisson random 
measures in this section. 

Theorem 4.2.1 There is a unique positive weak solution to (14.2.11) and the solution is a 
CBI-process with transition semigroup {Q])t>o- 

Proof. Suppose that {y{t)} is a cadlag realization of the CBI-process with transition 
semigroup given by (12.1.201) and (13.1.91 ). By Theorem 14.1. 21 the process has no negative 
jumps and the random measure 




(4.2.1) 



N{ds,dz) := ^l{y(^s)^y^s-)}S{s,y{s)-y{s-)){ds,dz) 



s>0 



has predictable compensator 




and 




(4.2.2) 



4.2. STOCHASTIC EQUATIONS OF CBI-PROCESSES 



49 



where N(ds, dz) = N(ds, dz) — N{ds, dz) and t H- M^{t) is a continuous local martin- 
gale with quadratic variation 2cy{t—)dt. By representation theorems for semimartingales, 
we have equation (14.2.11) on an extension of the original probability space; see, e.g., Ikeda 
and Watanabe (1989, p. 90 and p.93). That proves the existence of a weak solution to 
(14.2.11) . Conversely, if {y{t)} is a positive solution to (14.2.11) . one can use Ito's formula 
to see the process is a solution of the martingale problem associated with the generator 
L defined by (14.1.11) . By Theorem 14.2.11 we see {y{t)} is a CBI-process with transition 
semigroup {Q])t>o- That implies the weak uniqueness of the solution to (14.2.11) . □ 

Theorem 4.2.2 Suppose that m^dz) = az'^'^'dz for constants a > and 1 < a < 
2. Then the CBI-process with transition semigroup {Q])t>o is the unique positive weak 
solution of 

dy{t) = y/2cy{t)dB{t) + ^ay{t-)dzo{t) - by{t)dt + dz^(t), (4.2.3) 

where {B(t)} is a Brownian motion, {zQ{t)} is a one-sided a-stable process with Levy 
measure z~^~°'dz, {zi{t)} is an increasing Levy process defined by {f3, n), and {B{t)}, 
{zoit)} and {zi{t)} are independent of each other 

Proof. We assume cr > 0, for otherwise the proof is easier. Let us consider the CBI- 
process {y{t)} given by (|4.2.1I) with {iVo(ds, dz, du)} being a Poisson random measure 
on (0, oo)^ with intensity az'^'^'dsdzdu. We define the random measure {A^(ds, d^;)} 
on (0, ooY by 

NiiO,t]xB)= / / 1m.-)>o}1b ^—- )Noids,dz,du) 

Jo Jo Jo ^ y(^y\s-)^ 

pt poo 

+ / / / l{y{s-)=o}'^B{z)No{ds,dz,du). 
Jo Jo Jo 

It is easy to compute that {A^(ds, dz)} has predictable compensator 

+ I I 1{,(._)=o}1b(2)^ 

JO z 

Thus {A^(ds, d^;)} is a Poisson random measure with intensity z~^"°'dsdz; see, e.g., Ikeda 
and Watanabe (1989, p.93). Now define the Levy processes 

Zo{t) = / / zN{ds,dz) and zi{t) = (3t + / / zNi{ds,dz), 
Jo Jo Jo Jo 



Jo 

t POO 
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where N{ds, dz) = N{ds, dz) — N{ds, dz). It is easy to see that 

^/ayis-)dzois) = / / ^/ayi-s-) zN{ds,dz) 

zNo{ds, dz, du). 



JO JO 

t pco py{s-) 



JO JO 



Then we get (14.2.31) from (14.2.11) . Conversely, if {y{t)} is a solution of (14.2.31) . one can 
use Ito's formula to see that {y{t)} solves the martingale problem associated with the 
generator L defined by (14.1.11) with m{dz) = az^^^"dz. Then {y{t)} is a CBI-process 
with transition semigroup {Q])t>o and the solution of (14.2.31) is unique in law. □ 



Theorem 4.2.3 The pathwise uniqueness holds for positive solutions to (14.2.11) . 

Proof. For each integern > define a„ = exp{— n(n+l)/2}. Thena„ — decreasingly 
as n — )> oo and 




z dz = n, n > 1. 



Let X gn{x) be a positive continuous function supported by (a„, fln-i) so that 




gn{x)dx = 1 



and gn{x) < 2{nx) ^ for every x > 0. For n > let 

r\A ry 

fn{z) = dy I gn{x)dx, z eR. 
Jo JO 

It is easy to see that \fn{z)\ < 1 and 

< \zX{z) = \z\gn{\z\) < 2n-\ z e M. 

Moreover, we have fn{z) \z\ increasingly as n — oo. Suppose that {y{t) : t > 0} and 
{z{t) : t > 0} are both positive solutions of (14.2.11) . Let at = z{t) — y{t) for t > 0. From 
(14.2.11) we have 



at = aQ — h I as-ds + ^/2c / (a/ z{s) — \J y{s)^dB{s) 
Jo Jo 

pt POO pzi^s—) 

+ / zNQ{ds, dz, du). 

Jo Jo Jy{s-) 
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By this and Ito's formula, 

fnictt) = /n(«o) -b f'^{as)asds + c fni^s) {^/z(s) - >/y(s)Yds 
Jo Jo 

Pt POO 

+ / asl{a,>o}ds / [/n(«s + z) - /„,(«,) - z f!^{as)]m{dz) 

POO 

asl{a,<o}ds / [/„,(«« - z) - /„,(«,) + zfl{as)]m{dz) 
Jo 

+ martingale. (4.2.4) 

It is easy to see that |/„(a + x) — /„(a) | < |x| for any a, x G M. If ax > 0, we have 

\fn{a + x)- Ua) - xf'M\ < (2|ax|) A (n^VH. 
Taking the expectation in both sides of (14.2.41) gives 



P[/nM] < P[/„(«o)] + l&l / n\as\]ds + c ^[fl{as)\as\]ds 

Jo Jo 

POO 

ds j {{2zV[\as\]) ^{n-^z^)}m{dz). 
I Jo 

Then letting n — > oo we get 

P[|z(t) - ym < P[|z(0) - y(0)|] + |6| f nHs) - y{s)Wds. 



By this and Gronwall's inequality one can see the pathwise uniqueness holds for (|4.2.1I) . 

□ 

Theorem l4.2.3l was first proved in Dawson and Li (2006). By Theorems l4.2.1l and l4.2.3l 
there is a unique positive strong solution to (14.2.11) : see, e.g.. Situ (2005, p.76 and p. 104). 
The pathwise uniqueness of (14.2.31) was proved in Fu and Li (2010). 



4.3 Lamperti's transformations by time changes 

The results of Lamperti (1967b) assert that CB-processes are in one-to-one correspon- 
dence with spectrally positive Levy processes via simple random time changes. Caballero 
et al. (2009) recently gave proofs of those results using the approach of stochastic equa- 
tions; see also Helland (1978) and Silverstein (1968). Suppose that is a branching 
mechanism given by (12.1.131) . Let {x{t) : t > 0} be a CB-process with x{Q) = x > 
and with branching mechanism is given by (12.1.131) . Let {Yt : t > 0} be a spectrally 
positive Levy process such that Yq = x and 



\ogP exp{iX(Yt - Yr)} = {t-r)^{-i\), \eM,t>r>0. (4.3.1) 
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Let r = inf{s > : K, = 0} be its first liitting time at zero and let Zt = Yt/^r for t > 0. 
The proofs of tlie following two theorems were essentially adopted from Caballero et al. 
(2009). 

Theorem 4.3.1 For any t > let z{t) = x{K{t)), where 

Kit) = inf |n > : ^ x(s-)ds = x{s)ds > t|. (4.3.2) 
Then {z{t) : t > 0} is distributed identically on D[Q, oo) with {Zt : t > 0}. 

Proof. By the result of Theorem 14.2.11 we may assume {x{t)} solves the stochastic 
integral equation 

x{t) = x + I v^2cx(s-)d5(s) - I bx{s-)<ls 
Jo Jo 

pt POO px{s—) 

+ 11 I ^iVo(ds,d2,dM), (4.3.3) 



^0 ^0 



where {B(t)} is a Brownian motion, {A^o(ds, d^;, du)} is a Poisson random measures on 
(0, oo)^ with intensity dsm{dz)du and iV'o(ds, dz, du) = A''o(ds, d^;, du) — dsm{dz)du. 
It follows that 

z{t)=x+ / ^2cx{s-)dB{s) - / bx{s-)ds 
Jo Jo 

rn{t) poo /•x(s—) 

+ / / zNoids,dz,du) 

Jo Jo Jo 



'0 Jo JO ^ 

x+V2^W{t)-b z{s-)dK{s) 



t foo i'x{k{s)—) 

Jo Jo 







zNo{dK{s),dz,du), (4.3.4) 



where 



W{t) = / " y/x{s-)dB{s) = [ y/z{s-)dB{K{s)) 

Jo Jo 
is a continuous martingale. From (|4.3.2I) we have 

l{z{s-)>0}dK,{s) = 1{^(,_)>0}2:(S-)~MS. 

Let To = inf{t > : z{t) = 0}. Since zero is a trap for {z{t)}, we have 
t ft 

z{s-)dK{s) = / l{2(s-)>0}ds = t A Tq. 

Jo 
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Then has quadratic variation process {W)(t) = t Atq, so it is a Brownian motion 

stopped at tq. It is easy to extend {W{t)} to a Brownian motion with infinite time. Now 
define the random measure {A^(ds, d^)} on (0, oo)^ by 

iV((0,t] X (a,6]) = / / \{,^s-)>o}No{dK{s),dz,du), 



Ja JO 



where t > and 6 > a > 0. It is easy to compute that {A^((0, t] x {a,b]) : t > 0} has 
predictable compensator 

7V((0, t] X (a, 6]) = / m{a,b]z{s—)dK{s) = / m{a,b]l{s<To}ds. 
Jo Jo 

Then we can extend {A^(ds, dz)} is a Poisson random measure on (0, oo)^ with intensity 
dsm(dz); see, e.g., Ikeda and Watanabe (1989, p. 93). From (14.3.41) we conclude that 
{z(t)} is distributed on D[0, oo) identically with {Zt : t > 0}. □ 

Theorem 4.3.2 For any t > let Xt = Zg(t), where 

e{t) = inf {n > : ^ Z;}ds = Z;^ds > t|. (4.3.5) 
Then {Xt : t > 0} is distributed identically on D[0, oo) with {x(t) : t > 0}. 

Proof. By the Levy-Ito decomposition, up to an extension of the original probability 
space we may assume {5^t} is given by 

pt POO p1 

Yt = X + V2cW{t) -bt+ / / zMo{ds,dz,du), 

Jo Jo Jo 

where {W^(t)} is a Brownian motion, {Mo(ds, d^, du)} is a Poisson random measures on 
(0, oo)^ with intensity dsm{dz)du and Mo(ds, dz, du) = Mo(ds, d^, du) — dsm{dz)du. 
It follows that 



pt poo pi 

Xt = x + V2^W{e{t))-be{t)+ / / / zMo{de{s),dz,du). (4.3.6) 

Jo Jo Jo 

From (14.3.51) we have 



e{t)= [ Zeis)ds= [ Xsds. 
Jo Jo 



Then the continuous martingale {W{9{t))} has the representation 

W{e{t)) = [ v^d5(s), t > 
Jo 
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for another Brownian motion {B(t)}. Now we take an independent Poisson random 
measure {Mi(ds, d^, du)} on (0, oo)^ with intensity dsm{dz)du and define the random 
measure 

No{ds,dz,du) = l{u<x,_yMo{d9{s),dz,X;}du) + l{u>x^_}M^{ds,dz,du). 

It is easy to see that {Ai'o(ds, d^, du)} has deterministic compensator dsm{dz)du, so it is 
a Poisson random measures. From (14.3.61) we see that {Xt} is a weak solution of (14.3.31) . 
That gives the desired result. □ 



Chapter 5 

State-dependent immigration structures 



In this chapter we investigate the structures of state-dependent immigration associated 
with CB-processes. For simplicity, we only consider interactive immigration rates. The 
models are defined in terms of some stochastic integral equations generalizing (14.2.11) . 
We prove the existence and pathwise uniqueness of solutions to the stochastic integral 
equations. Similar immigration structures were studied in Li (2011) in the setting of 
superprocesses by considering different types of stochastic equations. We shall deal with 
processes with cadlag paths as in Li (201 1). 

5.1 Time-dependent immigration 

In this section, we introduce a generalization of the CBI-process. Let is a branch- 
ing mechanism given by (12.1.131) and let {Qt)t>o be the transition semigroup defined by 
(12.1.151) and (12.1.201) . We consider a set of functions {^s '■ s > 0} C given by 



where (3^ > and (1 A u)ns{du) is a finite measure on (0, oo). We assume s "ipsiz) is 
locally bounded and measurable on [0, oo) for each z > 0. By Theorems ll.2.3l and ll.2.4[ 
for any t > r > there is an infinitely divisible probability measure 'jr^t on [0, oo) defined 




z>0 



(5.1.1) 



by 




A > 0. 



(5.1.2) 



Then we can define the probability kernels {QJt ■ t > r > 0)hy 
Q^ti^, ■) := Qt-rix, ■) * 7r-,t(-), X > 0. 



(5.1.3) 
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It is easily seen that 

2 e-^yQl,{x,dy) = exp [ - xvt^riX) - M^t-s{X))dsy (5.1.4) 

Moreover, the kernels {Qj^t ■ t > r > 0) form a transition semigroup on R+. A Markov 
process with transition semigroup {QJt : t > r > 0) is called a special inhomogeneous 
CBI-process with branching mechanism cf) and time-dependent immigration mechanism 
{il)s : s > 0}. One can see that the time-space homogeneous transition semigroup associ- 
ated with {Qlt : t > r > 0) is a Feller semigroup. Then {Qlt '■ t > r > 0) has a cadlag 
realization X = 2/(t), Q^^.). In particular, if uns{du) is a locally bounded 

kernel from [0, oo) to (0, oo), one can derive from (12.2.41) and (15.1.41) that 

/•oo /*t 

yQUx,dy)=xe-'^'-''^ + j e'^^^-^V^Wds. (5.1.5) 

where 

POO 







The reader may refer to Li (2002) for the discussions of general inhomogeneous immi- 
gration processes in the setting of measure-valued processes. 



5.2 Predictable immigration rates 

Let be a branching mechanism given by (|2. 1.131) and an immigration mechanism 
given by (13.1.101) . In this section, we give a construction of CBI-processes with ran- 
dom immigration rates given by predictable processes. Suppose that (fi, P) is 
a filtered probability space satisfying the usual hypotheses. Let {B{t) : t > 0} be an 
(^t)-Brownian motion and let {poit) : ^ > 0} and {pi{t) : t > 0} be (^t)-Poisson 
point processes on (0, oo)^ with characteristic measures m(dz)du and n{dz)du, respec- 
tively. We assume that the white noise and the Poisson processes are independent of each 
other. Let A'^o(ds, d^, du) and A''i(ds, d^, du) denote the Poisson random measures on 
(0, oo)^ associated with {poit)} and {pi{t)}, respectively. Let iVo(ds, d^, du) denote the 
compensated measure of Ai'o(ds, d^, du). Suppose that p = {p{t) : t > 0} is a positive 
(^4) -predictable process such that t i-> P[p(t)] is locally bounded. We are interested in 
positive cadlag solutions of the stochastic equation 

ft poo /"Y,- 



Pt Pt poo r^a- 

Yt = Yo + a ^/Y^dB{s)+ / / zNo{ds,dz,du) 
Jo Jo Jo Jo 

pt pt POO pp{s) 

+ / {l3p{s)-hYs^)ds+ / / / zNi{ds,dzAu). (5.2.1) 
Jo Jo Jo Jo 
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Clearly, the above equation is a generalization of (14.2.11) . For any positive cadlag 
solution {Yt : t > 0} of (15.2.11) satisfying P[lo] < oo, one can use a standard stopping 
time argument to show that t P [Yt] is locally bounded and 

P[Yt] = P[Yo]+^'{0) [ P[p{s)]ds-b [ P[Y,]ds, (5.2.2) 

Jo Jo 

where tp'{0) is defined by (13.1.121) . By Ito's formula, it is easy to see that {Yt : t > 0} 
solves the following martingale problem: For every / G C^(M+), 

f{Yt) = f{Yo) + local mart. - b f f\Y,)YAs + ]-<j^ f f"{Y,)YAs 

Jo ^ Jo 

pt POO 

+ / YAs / [/(n + z)- f{Ys) - zf\Y,)]m{dz) 
Jo^ Jo 

+ ^ pis){pf'iYs) + J^^[fiYs + z)-fiYs)Hdz)]ds. (5.2.3) 

Proposition 5.2.1 Suppose that {Yt : t > 0} is a positive cadlag solution o/ (15.2. II) and 

{Zt : t > 0} is a positive cadlag solution of the equation with p = {p{t) : t > 0} 
replaced by t] = {ri{t) : t >0}. Then for any t>Qwe have 

P[\Zt -Yt\]< e'^'*{p[|Zo -Yo\]+ ^'(0) f^PMs) - p{s)\]ds], (5.2.4) 
where ifj'if)) is defined by (13.1.121) . 

Proof. For each integern > define a„ = exp{— n(n+l)/2}. Thena„ — )■ decreasingly 
as n -> oo and 

z~^dz = n, n > 1. 

71 

Let X (-)■ gn{x) be a positive continuous function supported by (a„, a„_i) so that 



gn{x)dx = 1 



and gn{x) < 2{nx) ^ for every x > 0. Let 



ry 

fn{z) = dy gnix)dx, z eR. 
Jo Jo 



It is easy to see that \fniz)\ < 1 and 

0<\z\f::{z) = \z\gM)<2n-\ z eR. 
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Moreover, we have fniz) — )• \z\ increasingly as n — t- oo. Let at = Zt ~ Yt for t > 0. 
From (15.2.11) we have 

at = ao + P [ [vis) - p{s)]ds + a [ {^ZT- - v^)d5(s) 

JO 

pt rt poo rZs- 

— b Os^ds + / / / zNQ{ds,dz,du) 
Jo Jo Jo Jys- 

mri(s) 
zNi{ds,dz,du). (5.2.5) 

By this and Ito's formula, 

fn{at) = /n(ao) + /3 / /n("s)[?7(s) - p{s)]ds -b /;^(a^)a^ds 
Jo Jo 

pt POO 

+ / «sl{a,>o}ds / [fn{as + z) - - z f^{as)]m{dz) 

Jo Jo 

- / asl{as<0}^S / [fn{0ts - z) - fn{ois) + Z f '^{u s)]m{dz) 

Jo Jo 

t 



+ / [v{s) - p{s)]l{^{s)>p{s)}ds / [fn{as + z) - fn{as)]n{dz) 



JO 

t POO 

[p{s) - r/(s)]l{p(,)>^(,)}ds / - z)- fn{as)]n{dz) 

Jo 
+ martingale. (5.2.6) 

It is easy to see that |/n(a + a^) — /n(a) I < l^^l for any a, x G M. If ax > 0, we have 

\fn{a + x)- Ua) - xf'M\ < {2\ax\) A {n^M^). 

Taking the expectation in both sides of (15.2.61) gives 

P[/nK)] <P[/„M] + /3 fv[\r,[s)-p{s)\]ds + \b\ ['p[\as\]ds 

Jo Jo 

z) + n^^aH 



pt poo 

+ / PMs) - p{s)\]ds / znid 
Jo Jo 

POO 

ds / {i2zP[\as\]) Ain^h'^)}m{dz). 
Jo 



By letting n — > oo we get 

P[\Zt-Yt\] < P[\Zo-Yo\] + \b\ fp[\Z,-Y.,\]ds 

Jo 
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+ ^'(0) / PUs) - pis)\]ds. (5.2.7) 
Then we get the desired estimate follows by Gronwall's inequality. □ 

Proposition 5.2.2 Suppose that {Yt : t > 0} is a positive cadlag solution o/ (15. 2. II) and 

{Zt : t > 0} is a positive cadlag solution of the equation with {h, p) replaced by {c.rf). 
Then for any t>Qwe have 



sup \Zs - Ys\ 

0<s<t 



<P[\Zo-Yo\]+^'{0) [ PUs) - p{s)\]ds 

t 

+ (|6|+2^ zm{dz)'^ J^P[\Zs-Ys\]ds 

+ 2a(^*P[|Z,-n|]ds)' 

+ 2 ( P[\Z, - Y,\]ds j^' z^m{dz)y, 

where ^'(0) is defined by (13.1.121) . 

Proof This follows by applying Doob's martingale inequality to (15.2.51) . □ 

Theorem 5.2.3 For any Yq > there is a pathwise unique positive cadlag solution {Yt : 
t > 0} of dSXTT) . 

Proof. The pathwise uniqueness of the solution follows by Proposition 15.2. 1 1 and Gron- 
wall's inequality. Without loss of generality, we may assume > is deterministic in 
proving the existence of the solution. We give the proof in two steps. 

Step 1. Let = ro < ri < r2 < ■ ■ • be an increasing sequence. For each i > I 
let Tji be a positive integrable random variable measurable with respect to ^r,_i- Let 
P = {p{^) : i > 0} be the positive (^j) -predictable step process given by 

oo 
i=l 

By Theorem 14.2. 1[ on each interval (rj_i,rj] there is a pathwise unique solution {Yt : 

ri-i <t <ri] to 

Yt = Yr^^,+a ^/Y^dB{s)+ / / zNo{ds,dz,du) 

ir,-! >''r,_i Jo Jo 
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+ / iPvi - bYs-)ds + / / / ^ zNi{ds,dz,du). 

Jri-i Jri-i Jo Jo 



Then {Yt : t > 0} is a solution to (ISXTl) . 

2. Suppose that p = {p(t) : t > 0} is general positive (^t) -predictable process 
such that t I— 7- P[p(t)] is locally bounded. Take a sequence of positive predictable step 
processes pk = {pkit) : t > 0} so that 



\pkis) - p{s)\ds 







(5.2.8) 



for every t > as k ^ oo. Let {Yk{t) : t > 0} be the solution to (15.2.11) with p = p^. By 
Proposition l5.2.1l Gronwall's inequality and (15.2.81) one sees 

sup P[\Yk{s)-Y,{s)\]^0 

0<s<t 

for every t>Oasi,A;— t-oo. Then Proposition 15 .2.21 implies 



sup \Yk{s) - Yi{s) 

0<s<t 







for every t > Oasi,/c — )• oo. Thus there is a subsequence {ki} C {k} and a cadlag 
process {F^ : t > 0} so that 

sup ^0 

0<s<t 

almost surely for every t > as i oo. It is routine to show that {Yt : t > 0} is a 
solution to (15.2.11) . □ 

Theorem 5.2.4 If p = {pit) : t > 0} is a deterministic locally bounded positive Borel 
function, the solution {Yt : t >Q} o/' (15. 2. II) is a special inhomogeneous CBI-process with 
branching mechanism and time-dependent immigration mechanisms {p{t)ip : t > 0}. 



Proof. By Theorem 14. 2. 11 when p{t) = p is a deterministic constant function, the process 
{Yt : t > 0} is a CBI-process with branching mechanism and immigration mechanisms 
pijj. If p = {p(t) : t > 0} is a general deterministic locally bounded positive Borel 
function, we can take each step function p^ = {pk{t) : t > 0} in the last proof to be 
deterministic. Then the solution {Ykit) : t > 0} of (15.2.11) with p = p^ is a special inho- 
mogeneous CBI-process with branching mechanism (p and time-dependent immigration 
mechanisms {pk{t)4' : t > 0}. In other words, for any A>0, t>r>0 and G G #r we 
have 



P[l^e~^n^'^] = P Icexpl -Yk{r)vt.r{X)- / Pfc(s)^K_.(A))ds 
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Letting k oo along the sequence {ki} mentioned in the last proof gives 



P[lGe 



-mi 



Icexp <^ - YrVt-ri\) 



Then {Yt : t > 0} is a CBl-process with immigration rate p = {p{t) : t > 0}. □ 

In view of the result of Theorem l5.2.4[ the solution {Fj : t > 0} to (15.2.11) can be called 
an inhomogeneous CBI-process with branching mechanism 0, immigration mechanism ip 
and predictable immigration rate p = {p{t) : t >0}. The results in this section are slight 
modifications of those in Li (201 1+), where some path-valued branching processes were 
introduced. 



5.3 Interactive immigration rates 

In this section, we give a construction of CBI-processes with interactive immigration 
rates. We shall use the set up of the second section. Suppose that z t-^ q{z) is a positive 
Lipschitz function on [0, oo). We consider the stochastic equation 

t ft POO pYg — 



Yt = Yo + a y/Ys-dB{s)+ / / zNo{ds,dz,du) 
Jo Jo Jo Jo 

ft ft foo fq{Ys-) 

+ [PqiYs.) - bYs.]ds + / / zNi{ds,dz,du). (5.3.1) 
Jo Jo Jo Jo 

This reduces to (14.2.11) when g is a constant function. We may interpret the solution 
{Yt : t > 0} of (15.3. 1|) as a CBI-process with interactive immigration rate given by the 
process s q{Ys_). 

Theorem 5.3.1 There is a pathwise unique solution {Yt : t > 0} o/ (15. 3. II) . 

Proof. Suppose that {Yt : t > 0} and {Zt : t > 0} are two solutions to this equation. Let 
K > Ohe a Lipschitz constant for the function z q{z). From (15.2.41) we get 

P[\Zt - Yt\] < ^^J'{0)e\'\' f P[\q{Zs) - q{Y,)\]ds < K^'{0)e\''\' f P[\Z, - Ys\]ds. 

Jo Jo 

Then the pathwise uniqueness for (15.3.11) follows by Gronwall's inequality. We next 
prove the existence of the solution using an approximating argument. Let lo(t) = 0. 
By Theorem l5.2.3l we can define inductively the sequence of processes {Yk{t) : t > 0}, 
k = 1,2, . . . as pathwise unique solutions of the stochastic equations 



Yk{t) =Yo-b f Yk{s~)ds + a f ^Yk{s-)dB{ 
Jo Jo 
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ft ft foo fYk{.s-) 

+ 13 q{Yu-i{s-))ds + / / zN^{ds,dzAu) 
Jo Jo Jo Jo 

ft foo fq{Yk^i{s-)) 

+ / zNi{ds,dz,du). (5.3.2) 



Jo Jo Jo 

Let Zk{t) = Yfc(t) - By (l5X4l) we have 



P[|Z,(t)|] < V^'(0)el^l* / P[|g(n_i(.))-g(n_2(5))|]d. 

Jo 

< i^^/.'(0)el''l* f P[\Zk-,{s)\]ds. 



By (15.3.21) one sees that {Zi{t) : t > 0} is a CBI-process with branching mechanism (p 
and immigration mechanism q{0)4'. In view of (13.1.1 II) . we have 

P[\Z^m = e-^*P[ro] + V^'(O) te-'^ds. 

Jo 

By a standard argument, one shows 

oo 

V sup P[\Yk{s) - < oo, 

so the Lipschitz property of z t-^ q{z) implies 

oo 

V sup P[|g(nGs)) < oo. 

It follows that 

lim fp[\q{Yk{s))-q{Yi{s))\]ds = Q. 

Then there exists a predictable process p = {p{s) : s > 0} so that 

lim / P[|g(n(s))-p(s)|]ds = 0. (5.3.3) 

Let {Yt : t > 0} be the positive cadlag process defined by (15.2.11) . By Proposition 15. 2.2[ 
there is a subsequence {/c„} C {k} so that a.s. 

lim sup |Yfc„(s) = 0, t>0. 

n^ooo<s<t 
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By the continuity of z q{z) we get a.s. 

limg(n„(s-)) = g(F(s-)), t > 0. 

n— >oo 

This and (|5.3.3I) imply that 

'V[|g(F(s-))-p(s)|]ds = 0. 



Then letting A; — )■ oo along {A;„ + 1} in (|5.3.2I) we see {F( : t > 0} is a solution of (|5.3.1I) . 

□ 

By Ito's formula, it is easy to see that the solution {Yt : t > 0} of (15.3.11) solves the 
following martingale problem: For every / G C^(]R_,_), 

fiYt) = fiYo) + local mart. - 6 / /'(r.)F.ds + ^-a^ [ f"iYs)Ysds 

Jo ^ Jo 

/*t /'OO 

+ / YAs / [/(n + z)- f{Ys) - zf\Y,)]m{dz) 
Jo^ Jo 

+ ^ + j^U^Ys + z)- /(n)]n(d^)}ds. (5.3.4) 

By Theorem 15. 3. 1[ the solution is a strong Markov process with generator given by 

1 /"°° 
Af{x) = -a\f\x) +xj^ [fix + z)~ fix) - zf'ix)]midz) 

-bxf'ix) + qix)[f3f'ix) + [fix + z)-fix)]nidz)y (5.3.5) 

We can also consider two Lipschitz functions z qiiz) and z h-> ^2(2) on [0, 00). 
By slightly modifying the arguments, one can show there is a pathwise unique solution to 

"t POO rYs- 



Yt = Yo + a ^/Y;:dBis)+ / / zNoids,dz,du) 
Jo Jo Jo Jo 

pt nt POO pq2{Ys-) 

+ [l3qiiYs^) - bYs^]ds + / / ziVi(ds,dz,du). (5.3.6) 
Jo Jo Jo Jo 



The solution of this equation can be understood as a CBI-process with interactive im- 
migration rates given by the processes s ^ qiiXs-) and s ^ q2(YsJ). This type of 
immigration structures were studied in Li (201 1) in the setting of superprocesses by con- 
sidering a different type of stochastic equations. 
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